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This article surveys a number of theoretical problems and open questions in the field of two- 
dimensional dilute Bose gases with weak repulsive interactions. In contrast to three dimensions, 
in two dimensions the formation of long-range order is prohibited by the Bogoliubov-Hohenberg 
theorem, and Bose-Einstein condensation is not expected to be realized. Nevertheless, the first 
experimental indications supporting the formation of a condensate in low dimensional systems 
have been recently obtained. This unexpected behaviour appears to be due to the non-uniformity 
introduced into a system by the external trapping potential. Theoretical predictions, made for 
homogeneous systems, require therefore careful reexamination. We survey a number of popular 
theoretical treatments of the dilute weakly interacting Bose gas and discuss their regions of appli- 
cability. The possibility of Bose-Einstein condensation in a two-dimensional gas, the validity of 
the perturbative t-matrix approximation and the diluteness condition are issues that we discuss 
in detail. 
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I. INTRODUCTION 

A. Revival of interest in low-dimensional systems 



Low-dimensional systems are interesting in general, as 
their low-temperature physics is governed by strong long- 
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range fluctuations. These fluctuations inhibit the forma- 
tion of the true long-range order (LRO), which is a key 
concept of phase transition theory in 3D. Thus, a 2D 
uniform system of interacting bosons does not undergo 
Bose-Einstein condensation at finite temperatures. How- 
ever, this system turns superfiuid below a certain tem- 
perature Tkt-i identified by Berezinskii, and Kosterlitz 
and Thouless (BKT) in 1971-73, signalling the presence 
of a so-called topological order. The elementary excita- 
tions of the superfiuid phase are pairs of vortices with 
opposite winding numbers. 

The experimental realization of such a system was 
for many years restricted to films of superfiuid 4 He 
on surfaces, which is also an example of a strongly- 
'interacting system. The breakthroughs in experimen- 
6 tal physics at the end of the last century have changed 
the situation drastically. The combination of laser cool- 
ing (S. Chu, C. Cohen-Tannoudji, W. D. Phillips, No- 
bel Prize for Physics, 1997) with evaporative cooling and 
magneto-optical traps provided experimental systems of 
cold atoms, which were primarily used to observe the 
long-awaited phenomenon of Bose-Einstein condensation 
(E. A. Cornell, W. Ketterle, E. Wieman, Nobel Prize for 
Physics 2001). The full tunability of magnetic and op- 
tical traps opens an extraordinary opportunity to study 
in practice not only ID and 2D Bose systems, but also 
dimensional crossovers under the influence of the num- 
ber of particles, size and shape of the system, interac- 
tion strength and temperature. These new developments 
have triggered a revival of theoretical interest in low- 
dimensional systems, when the old theoretical predictions 
are to be tested or carefully revised in order to address 
finite-size experimental systems, and a large field of new 
phenomena are to be explained. 

While the first experimental indications of the BKT 
transition in weakl y- interacting Bose system have been 
recently obtained ((Stock et all l2005(l . many questions 
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remain unanswered. One of the most interesting is, 
whether topological order survives under some conditions 
in the inhomogeneous trapped system, or is it dominated 
by the true LRO and Bose-Einstein condensation pre- 
vails? Can we control and directly observe the formation 
of vortex pairs in 2D quantum gases? These and other 
problems serve as the main motivation for this Collo- 
quium. 

In the next section we present a succinct overview of 
the history of work with dilute Bose systems, outlining 
some of the important theoretical problems relevant to 
weakly-interacting Bose gases. 



B. Historical overview 

The condensation of conserved particles that obey the 
same statistics as photons was predicted by Einstein in 
1924 even before the concept of Fermi statistics was in- 
troduced (1926). Einstein's prediction was preceded by 
an ingenious conjecture of Bose, who realized that black 
body radiation can be treated as a gas of indistinguish- 
able photons. Einstein generalized ideas of Bose to mate- 
rial particles and published two famous papers, in which 
he develo ped w hat we now call Bose-Einstein statistics 
!|EinsteinLll92lll92^ . 

The ideal gas of Bose particles is remarkably the only 
example of a non-interacting system in condensed matter 
physics that undergoes a phase transition upon decreas- 
ing the temperature. However, experimental realization 
of ideal Bose-Einstein condensates is extraordinarily dif- 
ficult, since realistic systems always involve interactions. 
Largely for this reason Einstein's ideas did not receive 
a wide recognition in the scientific community for many 
years as being devoid of any practical significance. The 
condensation phenomenon did not even appear in the 
textbooks, until in 1938 F. London recognized the anal- 
ogy between su perflu idity of liquid 4 He , disco vered by 
iKapitzal (^938) , and lAllen and Misenerl l)l938fl and an 
ideal Bose gas and emphasized that Einstein's statement 
was "erroneously discredited" l)Londorill938[) . 

In support of London's phcnomcnological ideas, the 
first microscopic theory of superfluidity in a system of 
weakly-interactin g Bose particles w as introduced in a 
brilliant paper bv lBogoliubovl l)l947j) . Subsequent discus- 
sion about the connection b etween superfluidity and Bose 
Einstein condensation led iPenrose and Onsagerl |)1956|) 
to formulate the generalized criterion for BE conden- 
sation. This line of research culminated in a paper of 

C. N. Yang, who in 1962 extended this criterion to su- 
perfluidity and superconductivity and proposed t he con- 
cept of off-diagonal long-range order (ODLRO) ( lYand . 
1962). The condensed phase is characterized then by a 
non-vanishing asymptotic of a one-body density matrix 
at large distances. 

During the decades which followed the work of Bo- 
goliubov, successful field-theoretical approaches were 
developed and many important predictions about the 



thermodynamics of the interacting Bose system were 
made. However, apart from the successful observa- 
tion of superfluidity in liquid Helium systems, the 
quest to create Bose-Einstein condensates (BEC) proved 
unrewarding for several decades. Finally, in 1995 
Bose-Einstein condensates were realized in a fascinat- 
ing series of experiments on rubid iu m and s odium 
vapours ^Cornell and Wiemari 120021 : iKetterleL l200ll 
iKetterle et all Il999j) 1. The importance of this experi- 
mental achievement was recognized in the 2001 Nobel 
Prize for Physics, shared by E. A. Cornell, W. Ketterle, 
and E. Wieman. 

The experimental realization of BEC has offered a 
unique opportunity to probe and control many inter- 
esting phenomena, not accessible or unstudied in the 
field of superfluidity, such as dimensional transitions, the 
crossover from Bose Einstein condensation to BCS pair 
condensation, interference effects, and disorder effects . 
Exotic links to c osmology llFedichev and Fischeit [2003). 
quantum optics ijRecati et al\ . l2005|) (two-state atomic 
quantum do ts within a condensate), and even w etting 
phenomena Ijlndekeu and Van Schaevbroeckl l2004|) have 
been recently proposed. The growing interest in Bose 
systems has resulted in more than 600 studies per year 
during last decade and the list of references related to 
BEC now exceeds 200 pages! 

The actual observation of condensation was hindered 
by enormous technical difficulties, so that even 15 years 
ago researchers dared not to believe that nature would 
ever provide them with the "right" system. The main 
problem to overcome is the condensation of most sys- 
tems into a solid or liquid upon cooling to low tempera- 
tures, which by-passes the BEC transition. In particular, 
the formation of clusters or molecules is driven by three- 
body collisions. The hard task for an experimentalist was 
therefore the creation of a gaseous system, in which three- 
body collisions occur much less frequently than two-body 
interactions. 

The gas in which the two-body interactions prevail is 
called dilute. Diluteness implies a very low density of 
the gas, so that the characteristic range a s of the poten- 
tial between the Bose particles is small compared to the 
mean particle distance, proportional to n" 1 / 3 in three 
dimensions (n — N/V being the density of the gas). The 
diluteness condition is therefore equivalent to the require- 
ment that the gas parameter be small 

n 1/3 a s < 1. (1.1) 

Ultra low density of the system leads to extremely low 
condensation temperatures (in the nanokelvin range), re- 
alization of which was another technical obstacle for the 
experimentalists. At low temperature the thermal ve- 
locity of the particles vt, which is proportional to the 
inverse De Broglie wave length 

/ 2irh 2 , , 

y mk B i 
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becomes very small (vt = h/mXr ~ 1 mm/sec) and at 
temperatures of the order of a few nK all the particles 
"jump" into a coherent ground state. Sufficient (Mute- 
ness of the gas is therefore one of the crucial conditions 
for BEC to be observed in the experiment. 

In order to reach the required regime of tempera- 
ture and density, vari ous cooling and trapping techniques 
have been developed l|Ketterle et alL \1999). Before being 
cooled atoms are confined in an external potential cre- 
ated by an applied magnetic field. The finite extent of 
the condensate cloud and its inherent inhomogeneity in- 
troduce a number of important differences between BEC 
in a trap and uniform gas. For example, a trapped gas of 
Bose atoms exhibits a BEC transition not only in momen- 
tum s pace, but in coordinate space as well ijDalfovo et all . 
1999). In practice however, condensates are so small that 
the literal observation of their size and shape is limited by 
the resolution of existing experimental equipment. Nev- 
ertheless real space Bose condensates provide a novel re- 
source for exploring many interesting phenomena, such 
as quantum interference effects and frequency dependent 
collective excitations. 

The effect of a magnetic trap becomes more dramatic 
for lower dimensionality of the system. For example, in 
2D a noninteracting trapped gas undergoes a BEC phase 
transition at finite temperature ifBagnato and Kleppnerl 
l!99lilLi et ~aH Il999l lwidomLll96Sj) in contrast to the 2D 
uniform case, where condensation is possible only at zero 
temperature. This difference arises due to modification 
of the density of states of the gas in the presence of a 
trap. 

The description of an interacting system in a 2D har- 
monic potential is not trivial. In the case of a uniform 
gas, long range order does not develop because of the 
preponderance of long wavelength phase fluctuations, in- 
herent to low-dimensional systems. This can be also seen 

as an infrared divergence of the integral J N(p) p^rfp , 
where N(p) is the number of particles out of the con- 
densate with momentum p. This divergence, on the 
other hand, is a consequence of the fact, that the energy 
of the system depends only on the phase gradient, and 
not on the phase itself, because the latter is n ot a well- 
defined quantity ijLifshitz and Pitaevskiill200^) . The ab- 
sence of long-range order in 2D systems with a contin- 
uous symmetry is often referred to as the Bogoliubov 
k~ 2 , or Hohe nberg-Merm i n- Wagner (BHM W) theorem 
(see works bv iBogoliubovl (1 1 96 ll Il99"l. IWagnerl ll 19661) . 
iMermin and Wagnerl ({1966(7 ! arid iHohenber 

i ©67)1, 

and we discuss this issue in m ore detail in Chapter llV.BI 
iFisher and Hohenberd l)l988j) pointed out that a conse- 
quence of the long-wavelength phase fluctuations is a 
drastic modification of the diluteness condition, so that 
the conventional low-density requirement for weakly- 
interacting 2D Bose gas, n 1 / 2 a s 1 is replaced by an 
inequality 

lnlnJ->l. (1.3) 

naj. 



Taken literally, condition (|1.3(l rules out the possibility 
of experimental realization of a 2D dilute Bose system. 
However, this condition does not work away from the 
transition. One can show from th e analy sis of quan- 
tum fluctuations (see iPetrov et all l)2004f) for review) 
that in this case the diluteness criter ion amounts to 
1/ ln(l /na 2 ) < 1, previously derived bv lSchickl l|l97lh . 

It is also intuitively clear that the trapping potential 
introduces a lower bound for the momentum of excita- 
tions and thus prevents the establishment of the long- 
range thermal fluctuations which destroy the con densate. 
Based on these arguments, IPetrov et all 1)200(1) showed 
the existence of a true condensate in a quasi-2D system 
in a wide parameter range. 

More generally, the BHMW approach is not really suit- 
able for a proper analysis of an inhomogeneous system, 
such as trap ped atomic v apour, a s pointed ou t bv lFische'rl 
l)2002ll2005() . In his work lFischeJ l)2002ll2005f) obtained a 
geometrical equivalent of the BHMW theorem, indepen- 
dently of the Hamiltonian of the system and showed that 
in the marginal d — 2 case true condensation is still pos- 
sible in an appropriately defined thermodynamic limit. 

In support of theoretical estimations, the first experi- 
mental confirmations of macroscopi c occupation of the 
harmonic oscillator gr ound state ijGorlitz et all . l200ll 
iRvchtarik et all . l2004|) became known in sodium atom 
vapours, confined to optical and magnetic traps. A rapid 
progress in experimental techniques made it possible to 
increase the aspect r atio (anis otropy) of the trap from 79 
llCorlitz et all . 1200 lh to 700 llSmith et all . l200ffl . This 
large anisotropy of the new traps is sufficient to con- 
fine condensates wi th ~ 10 5 atoms in a quasi-2D regime 
llSmith et flll . l2005D . Signs of local coherence were also 
observed in a two-dimensional g as of hydrogen atoms , 
absorbed on liquid 4 He surface l|Safonov et all Il998|) . 
Quasi-2D condensa te have been also recently created by 
IStock etal\ J2Q05) and interesting phase defects have 
been measured. The crossover from 3D condensates to 
two- and ultimately ID can be observed by changing the 
aspect ratio of the trap. 

As indicated in previous section, the recent progress 
in laser-based trapping techniques and creation of opti- 
cal lattices has led to a new generation of remarkable 
experiments. With controllable interparticle interaction 
it is now possible to observe the t ransition fro m the su- 
perfluid state to a Mott insulator (jBlochl l2004[) . Optical 
lattices provide a way to investigate various intriguing 
aspects of low-dimensional systems as well. Interest in 
2D configurations of Bose particles has arisen in the con- 
text of high-temperature superconductivity and the frac- 
tional quantum Hall effect. All in all, ultracold atomic 
gases have the potential to impact a very broad range of 
physics. 

In this Colloquium we discuss a number of selected 
issues related to two-dimensional weakly-interacting 
neutral Bose gases. If necessary, 3D problems 
are mentioned. We attempt to cover many refer- 
ences and otherwise refer the reader to numerous re- 
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sources, such as several excellent theoretical reviews 
llCastiill200HlDalfovo et a4ll999tlFetterl.l2002tlLeggettl 



l200lt iPetrov ^ all 12004 lYukalovl l2004l) 
iGriffinl 1199.-4 iPethick and Smithl l2002t I 



and books 
iPinesL [1962: 

iPitaev skii and Stri ngari l2003|h a Resource Letter for 
BEC (Hal| l2003h and BEC web-sites. Though a certain 
level of subjectivity is unavoidable, we aim to provide the 
necessary information about the field to those who feel 
lost after a preliminary contact with current literature 
but want to learn more about the main problems in the 
fascinating area of Bosc-Einstcin condensates in 2D. 



II. IDEAL BOSE GAS 

Consider a macroscopic system of noninteracting Bosc 
particles at finite temperature in the grand-canonical en- 
semble. The total number of particles in such a system 
is defined by the equation 



N 



E 

k 



p(e)n B (e)de, 



(2.1) 



where ns^k) — l/(exp ( S(e/ s — //) — 1) is the Bose-Einstein 
distribution function, /3 — l/fc^T and p(e) is the density 
of states. 

The chemical potential fi of the Bose gas, being nega- 
tive, increases as the temperature drops and vanishes at 
the critical temperature T c , indicating the phase transi- 
tion to a condensed state. The transition temperature is 
therefore defined by Ea. l|2.1|l with [i = 0. In d dimen- 
sions the density of states p(e) = dN e /de ~ e ( d - 2 )/ 2 5 and 
the particle density is proportional to the integral 



TV 
V 



S-^i^dt 

exp(e/T c ) - 1 ' 



(2.2) 



In 3D this integral converges and the Bose-Einstein con- 
densation temperature has a finite value, T 3D ~ n 2 / 3 . 
This result can be also understood as a temperature scale 
at which the thermal wave length becomes comparable 
with the average interparticle spacing At ~ I ~ n -1 / 3 . 
As A is proportional to T -1 / 2 , T c is proportional to n 2 / 3 . 

One can also calculate the and the number of particles 
occupying the ground state 



Y [ 1 - — 

T 



3/2 



(2.3) 



It is readily seen that No increases as the temperature 
decreases. This phenomenon of macroscopic occupation 
by particles of the state with minimal energy at low tem- 
peratures is referred to as Bose-Einstein condensation. 
Note that the actual condensation occurs in momentum 
space. 

In 2D, a constant density of states leads to an infrared 
divergent integral in expression 11'. 21 and condensation is 
not possible at any finite temperature. 



Let us discuss how this picture changes in the presence 
of a trap. The general treatment of this p roblem was con- 
sidered bv iBagnato and Kleopnerl l)l99l|) . They studied 
the possibility of the Bose-Einstein condensation of an 
ideal gas, confine d by one- or two-dimensional p ower-law 
trap: V ext ~ x v . IBagnato and Kleprmerl l|l99ll) showed 
that a two-dimensional system undergoes BEC for any 
finite value of rj, moreover, T 2D has a broad maximum 
in the vicinity of r/ = 2, i.e. for a trapping potential close 
to parabolic. (A one-dimensional system displays BEC 
only for 77 < 2.) 

Practically the confining trap is well approximated by 
a harmonic potential 



r ( , ( r» = |(^ 2 + ^v+^V). (2.4) 



For non-interacting particles we can write the many- 
body Hamiltonian as a sum of one-particle Hamiltonians 



H mb = Xa=i Hsp{i), whose eigenvalues are 
e« x ™ B n, = {n x +\)tw x +(n y +-)ftu) y +{n z +\)fw z . (2.5) 



The lowest energy of the system in the trap is e oo = f ftw, 
where for the sake of simplicity we introduced the average 
frequency ZU = (uj x + uj y + u) z )/3. 

Note, that in the ground state all N particles occupy 
the level eooo and the wave function of the "cloud" of 
these particles is easy to find 

0(ri..r N ) = JI<A>(rO 

i 

/muJho\ 3 / 4: ( m, 2 2 2\\ 

<A)(ri) = y — r— J exp[- — (u; x x + u; y y +cj z z)) 

(2.6) 



where 



VhO = (0JxU y U> z ) 1/3 . 



(2.7) 



In this case the density distribution of the particles is 
position dependent 



n(r) - N\tp (r)\' 



(2.8) 



and the first important length scale appearing in the 
problem is the size of the cloud 



O-hO 



m,Lj h o 



(2.9) 



which is just the average width of the Gaussian distribu- 
tion H2.6J1 (Fig.l). Experimentally aha is usually of order 
of lfim. 

Since we are mostly interested in low-dimensional ef- 
fects, it is instructive to mention the experimental re- 
alization of a two-dimensional atomic trap. An axially 
symmetric harmonic potential can be written in the form, 
V ext {r) = \m(jj\r\ + \mJlz 1 = ^muj']_(r']_+\ 2 z 2 ), where 
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FIG. 1 Column density of a cloud of trapped noninteracting 
bosons along the z-direction. The total density is a superpo- 
sition of a condensate density no and a thermal distribution 
of noncondensed particles ryr ■ 



X = lu z /lu± characterizes the degree of anisotropy. For 
ksT Huj z and fc^T > fko± the motion of atoms along 
the z direction is frozen (particles only undergo zero point 
oscillations), and kinematically the gas can be considered 
as two dimensional. Thus by making one dimension of 
the trap very narrow, oscillator states become widely sep- 
arated, and an effective 2D system is realized. 

At finite temperature only a fraction of the particles Nq 
occupies the lowest energy level and the others are ther- 
mally distributed over higher energy levels. However, we 
still can treat Nq as a macroscopic number. Thermal ex- 
citations will cause the size of the atomic cloud to grow 
with temperature. In the semiclas steal approximation 
ksT >> hujho, where the relevant excitation energies are 
much larger than the interlevel spacing, it can be shown 
that the size of the cloud increases as a square root of 

temperature Rt — a ho\J~^fj^- The important conclu- 
sion of this short discussion is that in harmonic traps, 
Bose condensation manifests itself as sharp peak in the 
central region of the density distribution in real space. 
The appearance of such a peak in both coordinate and 
momentum space is a peculiar feature of the trapped con- 
densates, with significant impact on both theory and ex- 
periment. This is very different from the uniform gas 
discussed above, where the condensation cannot be re- 
vealed in real space, for the condensate and uncondensed 
particles occupy the same volume. 

The total number of particles in the trap is defined by 



exp(- 



(2.10) 



which is derived from Ea. (|2.1|) with a discrete energy 
spectrum Q2.5|) . Note, that in this case the chemical po- 
tential at the transition point acquires a non-zero value 



of the lowest energy level: /i(T — > T c ) — > p c = (3/2)?SU. 

In the semiclassical approximation we can simplify 
(12.10(1 by replacing the summation with integration and 
a straightforward solution for fj, = fj, c gives the Bosc- 
condensation temperatures for the trapped gas in three 
and two dimensions 



T; 



3D 



(C(3))V3 



Who 



JVV3 



rp2D _ 



oN 1 / 2 



(2.11) 



(2.12) 



The two-dimensional condensation temperature is now 
finite (nonzero). This is related to the density-of-states 
effect of the gas in the trap. Indeed, in the semiclassi- 
cal approximation we can introduce a coordinate system 
defined by the three variables £ x ,y,z — ^>nx,y,z^>x,y,zi m 
terms of which the surface of constant energy l|2.5|l is the 
plane e = e x + e y + e z . Then the number of states A(e) 
is proportional to the volume in the first octant bounded 
by this plane 



A, 



1 



hO JO 



dr.. 



dc h 



de z = 



6hh 



hO 



(2-13) 

The density of states p — dN e /de is then quadratic in 
energy p 3D ~ e 2 in three dimensions and linear in energy 
in two dimensions p 2D ~ e, in contrast to the constant 
density of states of a uniform 2D gas, and the integral in 
the Ea. ((2.10(1 for /i = /i c is not infra-red divergent until 
d=l. 

It is now straightforward to calculate the condensate 
fraction (e.g. 3D) 



Aq 
A 



1 



/ T 



(2.14) 



and total energy of the system and correspondingly 
all the interesting therm odynamic quantities. In 2D 
the condensate frac tion is (jBagnato and Kleppnerl Il99ll 
iPetrov et all . 12004^ 



Ao 
A 



1 (2™) 



(2.15) 



Sign "w" in the expression 1(2.15(1 is related to the fact, 
that at T = T 2D the condensate fraction is not ex- 
actly zero, because there is a small correction to the 
resul t due to the finite number of particles in the sys- 
tem (jPetrov et all 12004^) . One should be therefore care- 
ful with the word "phase transition" in the context of 
trapped gases, because they are finite size systems and 
the phase transition notion is strictly defined only in ther- 
modynamic limit. It is better to say that at T c there 
is a sharp crossover to the BEC state in the system. 
Note also, that at T 2D the de-Broglie wavelength Ay be- 
comes comparable with the mean interparticle separation 
~ VT c /(Am^ ). 
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We end the section by remarking on the proper defi- 
nition of the thermodynamic limit in the trapped case. 
It is well known that the transition temperature should 
be well-defined in the thermodynamic limit. The usual 
definition when the ratio N/V is kept constant while the 
number of particles N and the volume V tend to infinity 
is apparently not suitable for the inhomogeneous situa- 
tion. The appropriately defined limit is then obtained by 
letting N — > oo and Uho —> 0, while keeping the prod- 
uct Nuj^q (or Nuif M in 2D) constant. In this case the 
temperatures (I2.11[) and 1)2. 12[) are well-defined. 

A comprehensive survey of various issues related to the 
behaviour of the ideal Bose gas in a harmo nic potential 
can be found in the paper bv lMullinl l|l997h . 

The ideal-gas results are summarized in the Table H] 

III. GROUND STATE OF A WEAKLY INTERACTING 
BOSE GAS 

A. Bogoliubov approximation 

In his seminal p aper "On the theory of superfluidity" 
l)Bogoliubovt Il947l) . published in 1947, Bogoliubov intro- 
duced the microscopic description of the ground state of 
a uniform, weakly interacting Bose gas. The assumption 
about the uniformity of the unperturbed ground state 
is crucial to his results. To assure a uniform Bose gas, 
Bogoliubov considered the case of repulsive interactions 
and made use of periodic boundary conditions. The gas 
is also assumed to be dilute (no 3 <S 1), which permits to 
simplify the many-body problem and account for inter- 
actions in a rather fundamental way. In contrast to the 
uniform case, the nonuniform ground state is very "sen- 
sitive" to the introduction of any interactions and makes 
the solution of the many body problem highly nontrivial. 

The standard Hamiltonian of an interacting Bose gas 

is 

H = i J v* + (0 V #(r)dr + 

+ i j \ + {r)t> + {r')U{r - r')^{r')^{r)drdr' (3.1) 

where U(r — r') is the interaction between particles. In 
momentum space this Hamiltonian reads 

H = ^2e p a+a p + U q a^a^,a p ^ q a p+qi (3.2) 

P pp'q 

U q = J e~ lqr U(r)dr is a Fourier component of 
the interaction, the bosonic field operator &(x) = 
l/y/V^2 p e tpx a p (here £ is a four vector), and the bo- 
son creation and annihilation operators satisfy the usual 
commutation relations [a p ,ai] = S pp i. 

Without interactions all N particles of the system oc- 
cupy the state with zero energy and zero momentum. 
The number of condensed particles iVo in this case is 
equal to the total number of particles N. When we switch 



on the interaction, two particles can scatter out of the 
condensate and occupy one of the the many zero-total- 
momentum states with separate momenta k and — k (in 
the lowest order perturbation theory) and Nq naturally 
decreases. 

For a dilute weakly interacting Bose gas one can as- 
sume that the total depletion of the condensate is small 
(SN/Nq <C 1) and most of the particles remain in the 
condensate Nq ^> 1. The key observation of Bogoliubov 
is that in this case the second-quantized condensate op- 
erators can be simply replaced by the "c"- number VNq 

d ,a^ ~ y/No. (3.3) 

The drawback of this prescription is that it leads to a 
Hamiltonian which no longer conserves the number of 
particles. This problem can be partly resolved by work- 
ing in the grand-canonical ensemble, in which additional 
terms —fiN p (N p — X) P /o a j^ a p) are introduced into the 
Hamiltonian (|4.27|) . This secures the conservation of par- 
ticles on the average. It is also worth mentioning that the 
Bogoliubov approximation is equivalent to the neglect of 
any dynamics in the condensed state. 

In the weak coupling limit the Hamiltonian 14.27fl can 
be diagonalized by applying the Bogoliubov canonical 
transformation 

a k = u k a k - v k at k 

al = u k a~l - v k a- k (3.4) 

and the resultant Hamiltonian describes the system of 
non-interacting quasiparticles with the spectrum 

/ ju2 £.4 

^ = y n oU - + — 2 , (3.5) 

where n Q = N /V is the density of condensed particles. 

From this dispersion relation 13.5fl it follows that in the 
long wavelength limit the Bogoliubov quasiparticles be- 
have as "phonons" with a sound velocity s — \J n °^° , and 
all of the low temperature thermodynamics of a Bose- 
condensed system is governed by this phonon spectrum. 
In the opposite limit of short wavelength the quasipar- 

.2 

tides behave as free particles with an energy j-. By 
equating the kinetic energy and the "Hartree" interaction 
energy uqUo one can straightforwardly find the "transi- 
tion" wave vector k c = ^/2mnoUo ~ \/2ms, which sepa- 
rates the phonon-likc behavior of elementary excitations 
from the free particle one. k c introduces an important 
length scale into the system (Fig. 2) 



A c = — = hj y/2mnoUo, (3.6) 

K c 

over which the coherence effects are important in the in- 
teraction between particles. It is usually called the heal- 
ing length (as in the context of trapped condensates), 
or sometimes the correlation or coherence length, and it 



7 



refers to correlations between excitations in the system. 
These correlations are distinct from the long-range cor- 
relations, which lead to condensation in the k = mode. 



phonons 
I 



crossover 



///////////// 




"free" particles 



k= JL 



FIG. 2 Length scale in the Bogoliubov problem: correlation 
length A c . k c ~ 1/A C separates the "free" particle behaviour 
from the linear dispersion region. Zigzag lines denote the 
residual interaction between the particles. 

One should also mention that the Bogoliubov canoni- 
cal transformation is equivalent to a summation over the 
most divergent terms in the perturbation-series expan- 
sion for the ground state energy. The summation of such 
series is also equivalent to making the random-phase ap- 
proximation (RPA). 

It was important in the theory of superfluidity, that 
the low-lying Bogoliubov quasiparticlcs follow a linear 
dispersion. This kind of behavior is fully consistent with 
the Landau criterion for superfluidity, i.e., that no exci- 
tation can be created in a liquid moving with a velocity 
v less than that of a sound (v < s). In case of noninter- 
acting particles the dispersion is quadratic for all k and 
superfluidity is not possible. 



B. Field-theoretical approaches: /-matrix approximation 

To go beyond the Bogoliubov approximation, one 
needs to take both multiple scattering diagrams and RPA 
contributions into account. That can be done for exampl e 
by means of a pseudopotential method ijLee et al\ . \l957\i . 
or by field-theoretical methods, fi rst appli e d to the Bose 
ga s of small density at T = b y lBeliaevI l|l958albl ) and 
bv lHueenholtz and Pines! |L959h 

The presence of the many particle condensate in the 
ground state of the interacting Bose gas was the main 
obstacle to application of the usual technique of Feyn- 
man diagrams to this system. Consider for example the 
one-particle Green's function in the interaction represen- 
tation 



G(x 



. <T{*(x)*+(x')S}> 

1 TW^ (37) 



(A*o = N). The S-matrix is expressed as usual 

S = T^exp^ — — J d^x\d^x 2 U (x\ — x%) x 

x ■^ + {x 1 )^ + {x 2 )^{x 2 )^{x l )^ (3.8) 

where X\ and x 2 are the four-vectors, and the interaction 
is U{x\ — x 2 ) = U{r\ — r 2 )5(ti — t 2 ). In order to de- 
rive the diagram series for the Green's function, we need 
to expand the S- matrix in powers of H int . Usually the 
terms containing the odd number of annihilation opera- 
tors vanish after averaging over the ground state, which 
unfortunately does not happen in the case of a Bose gas 
due to the above mentioned peculiarities of the ground 
state. The expectation value of the N product contain- 
ing ao apparently does not disappear and the standard 
method of constructing diagrams cannot be applied in 
the case of an interacting Bose gas. 

This difficulty was successfully resolved by Beliaev in 
f958. He noticed that for a large number of particles N 
the diagrammatic approach can be applied to particles 
with momenta p ^ 0, while the condensed phase (which 
does not disappear when the interactions are turned on) 
can be described as a sort of external field. It is thus 
convenient to separate the operators ao and aj (which 
act only on the ground state) from *F and , F + 



<F = <F' + a /W; <F+ = + a+ /VV. 



(3.9) 



The Green's function l|3.7|l is then also divided into two 
parts, and the operations T and < ... > are represented 
as two successive operations, the former acting only on 
^' and 'F + , and the latter acting only on ao and . 
The operators ao and aj, occurring in the S matrix are 
treated as parameters, and the expectation values over 
'F', 'F + ground state can be now calculated, using stan- 
dard techniques. 

With these ideas in mind, Beliaev succeeded in deriv- 
ing a general expression for the one-particle Green's func- 
tion of the interacting system in terms of some effective 
self-energies and the chemical potential /j,. However, 
the exact calculation of the Green's functions proved to 
be very complicated, and approximate methods of sum- 
ming the series of Feynman graphs were developed. 

For simplicity, Beliaev considered a short-range, cen- 
tral interaction potential U p = Uq for p < 1/a and 



Uu 



for p > 1/a. In the low density limit uqq <C 1, 



Here the average is taken over the ground state of N non- 
interacting Bose particles, which are all in the condensate 



where no is the density of the particles in the condensate, 
he obtained a crucial result, that the main contributions 
to the self-energies of the Green's function originate from 
ladder diagrams. In this case the real interaction U is 
replaced by an effective two-particle interaction T, rep- 
resenting the sum of contributions from all ladder type 
Feynman graphs (Fig. 3). The integral equation for the 
vertex T, called the Bethe-Salpeter equation is 
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T(x 1 ,x 2 ;x 3 ,x 4 ) = U Xl - X2 5{x 1 - x 3 )6(x 2 - x 4 ) 

+i / d 4 x 5 d 4 x 6 U xl ^ X2 G°(x 1 -x 5 )G°(x 2 ~ x 6 )T(x 5 , x 6 ; x 3 , x 4 ), 



(3.10) 



where x = (r, t) . In momentum representation this equa- 
tion reads 

r(pi,p2;?>3,P4) = U P1 - P2 + 

i J d 4 p 5 d 4 p 6 Up 1 -p 5 G (p 5 )G (p 6 )r(p 5 ,p 6 -,p 3 ,p i ), 

(3.11) 

where the momentum conservation condition p% + p 2 — 
P3 +P4=P5+P6 is implied and pt = (pi,p l )- 



Xi x 3 



x,=x 3 



>- 







>: > 






u 




+ U 


r 




>- 


> 


> > 




> 



x 2 x 4 



X 2 =X4 



X 2 x 6 x 4 



FIG. 3 Bethe-Salpeter equation for the two-particle scatter- 
ing vertex V. 

It is convenient to introduce relative and total mo- 
menta according to 



P1+P2 = P'\ P3 + Pi = P 

Pi-p 2 = 2p'; p 3 ~p 4 = 2p. 



(3.12) 



This transformation leads to the following equation 
t(p',p,P) = U(p'-p) + 

+l J C/ (P / -9) GO ( P / 2 + 9) G °( P / 2 -9M9 ) f^) 

(3.13) 

where we denote T in the center of mass representation 
by "£", and the free particle Green's function is G°(p) = 
(p° -p 2 /2m+ iSy 1 . 

A conventional t matrix equation is obtained from 
(|3.13l) after carrying out the integration over go. In two 
dimensions this results in the following equation 



t(p',p,P) = U p ^-J j^U p ,^ 



f (q, P, P) 



q — q 2 /m + iS 7 
(3.14) 

where fcg = P° — P 2 /4m. In the scattering theory this 
equation is also known as the Lippmann-Schwinger equa- 
tion. Physically the i-matrix corresponds to the renor- 
malization of the interaction by multiple scattering of one 
particle off another. 



The standard way to treat the dilute Bose gas is thus to 
replace the real potential, which is usually strongly sin- 
gular, by the zero momentum t matrix generated from 
multiple two-particle scattering, represented by the infi- 
nite summation of the ladder diagrams described above. 

The t matrix equation l|3.14|l cannot be solved explic- 
itly, but in general its solution can be expressed in terms 
of the scattering amplitude of two particles in vacuum. 
The scattering amplitude /(p',p) for a transition from 
the initial relative wave vector p to a finite relative vector 
p' is defined by an expression 



/(p',p) = J dqll(p' - q)* p (q), 



(3.15) 



where ^ p is a wave function of a scattering problem with 
potential U that satisfies the following Schrodinger equa- 
tion in momentum representation 



(k 2 



p 2 )^k(p) - J dqU(p - q)**(q) = 0. 



(3.16) 



According to elementary scattering theory CastinlEoOlt 
iDalfovo et all Il999t iFetterl . l2002t iLeggettl l200l|) . at low 



energies s-wave scattering becomes dominant, and the 
scattering amplitude fo is approximated to leading order 
by 



fo 



(3.17) 



(where the momentum dependence of the scattering am- 
plitude can be ignored in the low energy limit). Thus 
at low energies, in vacuum the only remaining parame- 
ter characterizing the interaction is the s-wave scattering 
length ag- 
in general, the t-matrix (|3.14|l requires the knowledge 
of the scattering amplitude for fcg 7^ 9 2 / to j known as 
"off-the-energy-shell" t-matrix. For two-particle scatter- 
ing in vacuum, discussed above, only on-shell t-matrix 
is physically relevant. In the situation when three-body 
collisions become importan t, the calculat ion of the off- 
shell t-matrix is necessary llFadeevl 1196(1) . In the con- 
text of the dilute Bose gases the off-shell t-matrix arises 
in conn ection with so-cal l ed many-body t-matr i x ap- 



proach llBiilsma and Stool [l997t iProukakis et all Il99: 
IStoof and Biilsmal . Il993|) . which we discuss in the next 
Chapter. The many-body t-matrix takes into account 
the effect of the medium (mean field) in which the col- 
lisions occur. At the low energy limit the many-body 
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t-matri x is approximated by the off-shell two-body t- 
matrix ijMorean et all . I2002T) . The solution of the off- 
shell t-matrix was first proposed bv lBeliaevI (|l958b^ and 
iGalitskiil The alternative approach based on 

the inhomogeneous Schrodinger equation, which allows 
to treat the hard-sphere central potentia ls in one, two 
and three dimensions , was considered by iMorean et all 
(2002). M organ et all l)2002T) have shown for any dimen- 
sion that for all potentials with a finite range, the long- 
wavelength limit of the off-shell t-matrix depends only 
on energy and not on the initial and final relative mo- 
menta of the scattered particles. This result means that 
low-energy collisions can be represented by a contact po- 
tential. 

Consider now the quasiparticle spectrum within the 
first order Beliaev approach. It turns out one can repro- 
duce the Bogoliubov result (|3.5[) with the only difference 
that instead of potential Uq the momentum independent 
scattering amplitude /o appears, for in the first order 
Uq is equal to Jq. The healing length i|3.6|) can then be 
related to a scattering length 



1 



TTa s n 



(3.18) 



The second order approximation does not modify the 
physical picture of the low temperature behaviour of the 
interacting Bose gas, but provides the corrections to the 
sound velocity, and a damping proportional to p 5 related 
to the process of decay of one phonon into two. The third 
order corrections involve the solution of a three particle 
problem, which to date has not been solved. 

We now turn to the 2 D system. Fo llowing the meth- 
ods developed by Beliaev. TSchickl l)l97l|) examined a two- 
dimensional system of hard-disk bosons of diameter a 
at low densities and absolu te zero (see also the recent 
study of lOvchinnikovl l)l993ft . The dimensionless expan- 
sion parameters are the interaction Uq and the gaseous 
parameter no 2 , which is small in the dilute limit. The 
application of Bcliacv's method to 2D systems is not as 
straightforward as it is for 3D systems. In the 3D case, 
the ladder diagrams are the only contributions which do 
not depend on the small parameter na 3 and therefore 
it is natural to take them into account while calculat- 
ing the first term in the expansions of all quantities in 
terms of density. In 2D the contributions from the lad- 
der diagrams depend logarithmically on the parameter 
na 2 , in particular, the effective interaction, or t matrix is 
proportional to 1/ ln(l /na 2 ) 



Jo 



2D 



4n 



m\n(l/na 2 



(3.19) 



The key conclusion of ISchickl l)l97l|) is that 
l/ln(l/na 2 ) plays the role of the small parameter in the 
2D dilute system at zero temperature and the dominant 
contributions are derived from the diagrams of first order 
in this parameter. In this approximation he calculated 



the leading order correction to the chemical potential 
4:Trh 2 n 



|l + 0[l/ln(na 2 )]} (3.20) 



mln(na 2 ) 

and the quasi-particle excitation spectrum 



k 2 fc 4 fc 4 

6c = \ u— + 7— = \ hrr 



Ann 



m 4m V 4m mln(l/na 2 ) 



k 2 . (3.21) 



In the long-wavelength limit the quasi-particles be- 
have as phonons with a speed of propagation s = 
\f— Airn/{m ln(na 2 )). The spectrum changes from 
phonon-like to free particle-like in the vicinity of the mo- 
mentum k c defined as 

ka <C k c a = — 167ma 2 (m ln(na )) <C 1 . 

The ground state energy per particle a nd the conden- 
sate fraction take the form ijSchickl Il97lfl 



E/N = - 



HQ 



m ln(na 2 ) 
1 

ln(na 2 ) 



{l + 0[l/ln(na 2 )]} 

+ [l/(ln(na 2 )) 2 l . (3.22) 



C. Gross-Pitaevskii mean-field theory 

The ground state and thermodynamic properties of an 
interacting Bose system confined to an external potential 
V ex t = \^m{r / aho) 2 (a/io is the trap size (|S30 can be 
directly calculated from the Hamiltonian 



H = J drijj+{r) (-— v 2 +V ext (r)j ^(r) + 

+ drdr'iP + {r)iP+{r')U{r - r')i/)(r')ip(r) (3.23) 

using numerical methods, such as quantum Monte Carlo. 
Nevertheless for most experimentally relevant situations 
(when the number of atoms is large) the mean-field de- 
scription of the system proves to be sufficient. In this 
case the macroscopic low energy behaviour of the sys- 
tem can be explored under the assumption that the or- 
der parameter varies over distances larger than the mean 
interparticle spacing. 

Such a mean-field approximation was first developed 
by Gross and Pitaevskii. Their approach, which is valid 
in the dilute limit, is a straightforward generalization of 
Bogoliubov theory for the gas in the trap. One should 
bear in mind that the diluteness condition n max a 3 <C 1 
does not automatically secure the weakness of the inter- 
actions. The interaction strength is specifi ed by an extra 
param eter (see in partic ular the revie w of lDalfovo et all 
(|1999|) and the paper bv lFetterl l)l999|n . The interaction 
energy, which is of the order of gNn is to be compared 
with the kinetic energy, proportional to Na^ 2 . Since 
the average density of atoms n ~ N/a^Q, the interaction 
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strength can be characterized by a dimensionless param- 
eter N\a s \/a,ho- When Na s <C a^o, it means that the 
coherence length A c (|3.6|) is large in comparison with the 
size of the trap aho and the system is assumed to be 
nearly ideal and is described by a Gaussian distribution 
()2.6(l . In the opposite limit, Na s ^S> a,ho, the coherence 
length is small a nd the dilute gas exhi bits important non- 
ideal behaviour l|Dalfovo et all Il999[) . 

The mean-field Gross-Pitaevskii approximation is ex- 
tensively presented in the lite rature (see for instance 
review bv iDalfovo et all (|l999j) . and paper by iLeeeettl 
(2003 ), and a review with an emphasis on experiment by 
lAngilellaeit al\ l)200fij) ). therefore we only mention briefly 
the key concepts of its derivation. Gross and Pitaevskii's 
approach is based on the Bogoliubov prescription for the 
condensate Ij3.3|l , according to which the boson field oper- 
ators if) are written as a sum of a classical field <f>, having 
the meaning of the order parameter, and a small pertur- 
bation if>' 



(3.24) 



implying that the depletion of the condensate is small. 
As a side note, we mention that in principle the problem 
of the order parameter definition in a finite inhomoge- 
neous system arises in this case, but it turns out that 
the wave function of the condensate has a clear mean- 
ing, if determined through the diagonalization of the one- 
body density matrix in analogy with liquid-helium drops 
l|Lewart et all . 1 19881). This issue is also discussed in de- 
tail in a review bv lLeggettl <|200lh . 

One can expand the theory in the parameter ip' 
and derive the equation for <fr either from the stan- 
dard Heisenberg equation, or alternatively by taking 
the variation of the classical action S of the type S = 



J dtdrcj) 



id t 



2m 



> with respect to </> 
The derivation of the 



(saddle point approximation) 
Gross-Pitaevskii (GP) equation and the next order cor- 
rections withi n the bosonic fi eld theory can be found in 
the paper bv lStenhoiml l)l998ft . 

The resulting Gross-Pitaevskii equation is 



2m 



V 2 +V ext +g\cj ) (r,t)\ z <t>(r,t), 



(3.25) 

where we have approximated the potential by a 5- 
function, V{r — r') = g5(r — r') (which we can do under 
the assumption that the interparticle spacing is much 

A ft 2 

larger than the interaction range), and where g = „ t as 
is the 3D coupling constant. This coupling constant is 
equal to the zero momentum limit of the scattering am- 
plitude H3.17fl discussed above. 

In the limit N — ► oo (Thomas-Fermi approximation) 
the kinetic energy contribution can be neglected and 
the Gross-Pitaevskii equation can be solved analytically. 
This classical Thomas Fermi approximation breaks down 
in the vicinity of the boundary of the condensate, where 
the gradient of the condensate density is no longer small. 



We discuss now the coupling constant of the 2D Bose 
gas. It was first demonstrated by Y . Lozovik in 1971 
(see the review of lPetrov et all l|2004|) ) that to zero or- 
der in perturbation theory the coupling constant g 2D = 
h 2 fQ D /m, where f§ D is the scattering amplitude at en- 
ergy of the relative motion E = 2\x. 

This coupling cons tant can be treated as a par ame- 
ter, as in the work bv l Bavindir and Tanatarl <|1998[l (see 
also references therein), where the two-dimensional Bose 
gases described by the GP equation have been studied. 
For some range of interaction strength it was shown that 
interacting bosons behave similarly to the noninteract- 
ing case in a harmonic trap. For weak short-range in- 
terparticle interactions, a finite temperature BEC phase 
transition was found to occur. 

On the other hand, the coupling constant in 2D is 
expected to display a logarithmic dependence o n den- 
sity (c f. (|3.19Jl ) in accordance with estimations bv ISchickl 
l|l97l|) for fo in case of a homogeneous gas. The precise 
choice of q 2D ha s in fact been a controversial issue (see 



Kim et 



2001) and references therein). For example, 
|1999j) suggested g 2D ~ l/ln(l/fca), where 
< ka <C 1 and k is the infrared cut off introduced by 
the trap at 1/ano, so that g 2D ~ 1/ ln(a/ l0 /a). This kind 
of approximation may be reasonable when the size of the 
trap is much larger than all other length scales in the 
problem. 

Note, that for qua si-2D gas in a trap th e coupling con- 
stant was derived bv lPetrov et all l)2000j) 



~Q2D 



2V2^h 2 



a ho /a + (l/V^) ln(l/TTk 2 a 2 h0 ) 



(3.26) 



The rigorous derivation of the Gross-Pitaevskii func- 
tio nal for a two-dime nsional interacting gas was provided 
bv lLieb et all l)200l[) . Their analysis leads to the follow- 
ing expression for the coupling constant 



,2£> 



\ln(na 2 )\ ' 



(3.27) 



where n is the average density of the particles, propor- 
tional to v^V- The mean density is defined as n — 
1/N J n TF (r) 2 d 2 r, with the Thomas-Fermi density be- 
ing n TF (r) = (/j. tf — V ex t{r))/8n, and /z TF chosen so 



that the constraint J 



-.TF 



N holds. The density ex- 



pansion has been appli ed to the case of a 2D Bose gas 
at zero temperature bv lChernv and Shanenko! l)200lj) in 
order to derive the Gross-Pitaevskii equation. 

The modification of the GP equation due to the many- 
body renormalization of the scattering , mentioned in th e 
section IIII.BI has been provided by iLee et all l)2002l) . 
The effective interparticle interaction in 2D is modeled 
by the off-shell two-body t-matrix, that at low energies 
depends on the energy of the collision. The energy de- 
pendence of the effective interaction can be written in the 
density-dependent form and applied to trapped 2D gas. 
This leads to the GP equation, describing the condensate 
wave-function that no longer has a cub ic non-linearity i n 
but instead goes as (|-0| 2 / ln IV'I 2 )^ <|Lee et all . 1200^1 . 
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It is also interesting to analyze the deviations from 
the mean-field behaviour, since the experimental system 
is well controlled nowadays and different regimes can 
be realized. The corrections to the mean-field ground 
state solution stem from the quantum fluctuations, and 
their effect becomes more prominent with the growth 
of the gas parameter, as has been observed in Monte 
Carlo simulations. For the calculation of quantum cor- 
rections i n a systematic way we r e fer th e reader to the 
paper bv I Andersen and Haueerudl l)2002|) and references 
therein. Many references on the GP approximati o n and 
beyond can be found e lsewhere ijAngilella et all l2004t 
iKolomeiskv et ad 120001. For the effects of a third spatial 
dimension and the self-consist ent calculation of the cou- 
pling constant see the paper of lChernv and Brandl ( 20041 
and references therein. 



IV. FINITE TEMPERATURE PROBLEMS 

Zero-temperature techniques are not really suitable for 
controlling IR thermal fluctuations, and new methods 
have to be devised to describe the interacting system at 
finite T. At the beginning of the Chapter HVI we present 
the generic properties of the 2D XY models and the con- 
cept of quasi-long-range order, which is the central con- 
cept in the phase transition theory in 2D. Why the true 
long-range order cannot form in 2D uniform system, we 
discuss in detail in section HV.BI and especially the way 
familiar concepts from 2D phase transition theory should 
be revised in the trapped case. 

In Section IIV.CI we present the theory of Popov, who 
pioneered the finite-temperature generalization of Beli- 
aev's field-theoretic approach, and described the low- 
temperature superfluid state of the 2D Bose gas. In Sec- 
tion llV.Dl we show how the diluteness condition of Fisher 
and Hohenberg, discussed in the introduction, arises as 
an applicability limit of the Popov's t-matrix approach. 
In Section TlV-EI we describe methods, which generalized 
and/or improve the results of Popov, and also the Monte 
Carlo simulations, which are up to date most reliable 
numerical calculations of the superfluid phase of the 2D 
Bose system. Before concluding, we mention how unique 
the 2D system with a contact interaction is, for it pos- 
sesses an inherent symmetry, which leads to the birth of 
the special breathing modes, which in principle can be 
checked experimentally. 



A. Introduction: 2D XY models 

For our further analysis it is important to recognize 
that a uniform, interacting Bose system belongs to the 
XY universality class, characterized by a vector order 
parameter (for a comprehensi ve analysis see the book by 
IChaikin and Lubenskvl l)l995l1 1. It means that the finite- 
temperature behaviour of the 2D Bose gas is determined 
by generic properties of the 2D XY model. 



We know, that the 2D XY models are very special, 
for the long-range thermal fluctuations destroy the long- 
range order at finite temperatures (Bose-Einstein con- 
densation in case of 2D Bose gas). The existence of 
these long-wavelength modes in a 2D Bose fluid was 
first pointed out by Bogoliubov in his "fc~ 2 theo rem" 
in 1961, and later confi r med by iHohenbergl l|l967j) and 
bv lMermin and Wa gner (1966J) (this issue is discussed in 
more detail in section llV. Bp . 

However, a special type of order - topological 
order - which gives rise to superfluidity, can de- 
velop in a two-dimensional Bose fluid below the 
"Kosterlitz-Thouless" temperat ure T kt, as predi c ted by 
Kosterlitz and Thoulessl <|1973|) and Berezinskiil (jl970t 
19711) using the renormalization group method (RG). Be- 
low Tkt the continuous U(l) symmetry (rotations in 
a two-dimensional plane) is broken and the system ac- 
quires a finite rigidity, or phase stiffness p s . The or- 
der parameter correlations decay algebraically (for any 
coupling of the XY model), and the average order pa- 
rameter is zero. However locally, the order parameter 
can have a well-defined value. This unique situation 
is described in terms of quasi- long-range order (QLRO) 
ijChaikin and Lubenskvl Il995j) . Important low- lying ex- 
citations of the QLRO phase are vortex pairs (two vor- 
tices with opposite winding numbers) whose fugacity de- 
creases with distance, thus not destroying the connectiv- 
ity of the state (therefore p s ^ 0). 

The phase transition to a disordered state (with p s = 
0) is associated with a dissociation of the coupled vortex 
pairs. Above Tkt the vortex fluid can be treated as a 
kind of vortex plasma, where vortices play the role of 
mobile "charges" , interacting via a Coulomb potential. 
In this language the state below Tkt can be described as 
an "insulating" state of bound "charges" . The mapping 
of the 2D XY model onto the two-dimension al Coulomb 
gas is considered in detail in the review by iMinnhaeenl 

The rigidity, or superfluid density p s does not go con- 
tinuously to zero at the critical temperature, but experi- 
ences a universal jump 

m 2 k B T K T _ vr 

h 2 p s (T KT _) ~ 2' [ ■ > 

first predicted bv lNelson and Kosterlitzl l)l977j) and suc- 
cessfully verified in exper iments on superfluid 4 He films, 
absorbed on a substrate ijBishop and Repp"vl . ll978|) . 

An interesting interpretation of Kosterlitz-Thouless 
physics in the contex t of b osonic systems was put for- 
ward bv lKagan et al\ <)l987j) almost 20 years ago. They 
propose that below Tkt the system forms a "quasi- 
condensate" , a condensed state achieved in a local sense. 
The introduction of the quasi-condensate concept was 
motivated by a peculiar behaviour of the one-particle 
density matrix p(r) at large distances in 2D (Fig. 4). 

There are two length scales associated with the be- 
haviour of p{r): the aforementioned correlation length 
A c at which p(r) relaxes from the value n at r = 
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FIG. 4 One-particle density matrix p(r) = (ip + (0)ip(r)} in 
two dimensions. Two characteristic length scales Rg and A c 
are shown, Rg S> A c . At large distances A c <C L <C Re the 
one-particle density matrix is equal to the condensate density 
no- 



to no, and the characteristic radius of the phase fluc- 
tuations Rg, which is rather large Rg 3> A c . The ap- 
pearance of large Rg can be understood in the follow- 
ing w ay: at large distanc e s p fa lls off as a power law 
of r <|Kane and Kadanoffl ^1967^ ^1 p(r) ~ n (r/r*)- a , 
where r 3> r* and the coefficient a is proportional 
to the temperature and the Schick's parameter a ~ 
T/(T*ln{l/na 2 )) (T* ~ (2mA;?)- 1 ) and therefore is very 
small a <C 1 • As a result of this the densit y matrix p de- 
cays over a large length scale Rg ~ r*e 1 / Q ijKagan et all 
1987). 

Conceptually, the system can be divided into blocks 
of size L, which is smaller than Rg. In each block 
one can introduce the wave-function of the condensate 
with a well-defined phase. The whole system is then de- 
scribed in terms of an ensemble of wave-functions of the 
blocks. Condensate wave-functions within the ensemble 
corresponding to blocks separated by a distance greater 
than Rg have uncorrclated phases, and it is impossible 
to define the condensate wave function for the entire sys- 
tem as a whole. The state of matt er with a fluctuatin g 
phase is called a "quasi-condensate" l)Kaean et alJ . ll98'r]) . 
See also the ex tention of Bogoliubov m ethods to quasi- 
condensates bv lMora and Castinl l)2003[) . 

What happens to the XY universality class concepts 
in an experimentally realizable system of cold atoms con- 
fined in a trap, remains a controversial issue. We will see 
in the next Section, that many issues should be crucially 
reformulated in order to address the physics of trapped 
cold gases. 



B. Problems of the long-range order formation in 2D 

The notion that the development of the long-range or- 
de r (LRO ) is no t possible in 2D dates back to the work 
of iPeierlsl l)l935|) . who argued that the thermal motion 
of low energy phonons will ruin the LRO in a 2D solid. 
A rigor ous proof of the Peierls' statement was provided 
later bv lMerminl l)l968|h 



Subsequent work bv iMermin and Waenerl (^966) pro- 
vided a proof that there is no spontaneous magnetization 
or sublattice magnetization in an isotropic Heisenberg 
model with finite range interactions. At the same time 
Hoh enberd l)l967|) succeeded in ruling out the existence 
of a conventional superfluid or superconducting order- 
ing in one and two dimensions. It was also shown by 
IColemanl (|l9T2l) , that "there are no Goldstone bosons in 
2D" , which is equivalent to saying that there is no LRO 
in 2D. 

A rigorous proof of the Mermin-Wagner-Hohenberg re- 
sults exploits the Bogoliubov and Schwartz inequalities 
(Appendix I) and leads to the following result for the 
average occupation number of k states 



{ala k ) = n k > -- 



1 mTno 



k 2 7 



(4.2) 



Here no is a condensate density and n is a total density. 
It is clear now that the appearance of the condensate 
(macroscopic occupation of a single state) in 2D for finite 
temperatures fails due to the mathematical fact that the 
function k~ 2 is not integrable at small momenta in two- 
dimensional k space. Physically, the long-range thermal 
fluctuations prevent the formation of a coherent conden- 
sate. 

The same result can be obtained from the infrared 
asymptote of the one-particle Green's function at zero 
frequency 



G(k,0) 



no to 
n s k 2 



(4.3) 



and was first derived bv iBogoliubovl l(l96ll) . The deriva- 
tion of the asymptotic behaviour (14.311 in a fun ctional 
integrals approach can be found in lPopovl l)l983|h bmce 
the Green's function defines the average number of parti- 
cles with momentum k, it is readily seen we arrive to the 
same result Q4.2p . The statement, that the condensate 
does not appear in a 2D interacting Bose system at any 
finite temperature, is also known as the Bogoliubov k~ 2 
theorem. 

We have already mentioned that in the context of mod- 
ern condensed matter theory the absence of the LRO 
in 2D is discussed in terms of general properties of the 
XY models. A respective direction or a phase of the d- 
dimensional XY order parameter is specified by an angle 
9. The variance in the fluctuation of the order parameter 
phase is given by the integral 



Ps 



d d q _ 2 
(2n) dq 



TA 



d-2 



p s {d-2Y 



(4.4) 



where A is th e w ave number cutoff 
l)Chaikin and Lubenskvl Il995tk It is readily seen 
that d = 2 is the critical dimension of the XY universal- 
ity class and the fluctuations destroy long-range order 
in the 2D XY model in accordance with the conclusions 
of Bogoliubov, Mermin, Wagner and Hohenberg. Quasi- 
long-range order, discussed in the previous section is 
nevertheless possible in 2D. 
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In case of a trapped gas the Bogoliubov-Mermin- 
Wagner-Hohenberg (BMWH) the orem ru l es out BEC in 
2D in the interacting system fsee iMuliml l|l997f )1. How- 
ever, the question arises if one can actually apply BMHW 
theorem to a system confined within a harmonic poten- 
tial. Is it still possible to unambiguously rule out the con- 
densate formation in 2D atom traps? The applicability of 
the BMWH theorem to the inhomogeneous case requires 
careful consideration, for the Bogoliubov-Hohenberg in- 
equality was derived assuming an infinite uniform sys- 
tem. In this approximation, many special features of 
practically realized condensates, such as their formation 
in real space, are excluded. 

An alternative version of the Hohenberg inequality, 
suitable for the experimental ly realizable Bose sy stems, 
has been recently proposed bv lFischerl (|2002. 2005). Tak- 
ing the dimension of the trap to be an experimentally 
controlled parameter, Fischer addressed the issue of a 
spatially localized Bose condensate, with the question in 
mind of how far one could "stretch" the 3D condensate 
cloud before the coherence will be destroyed. Fischer de- 
rived an inequality, which controls the size of the smallest 
possible condensate for a given condensate and density 
profile. In Appendix I we briefly sketch the underlying 
concepts of his derivation. 

The resulting inequality reads 

> C(k) - JL (1 - |*(k)|^o)) , (4.5) 

where ifio(k) is the condensate wave- function, R c is the 
effective radius of the condensate wave function (effective 
radius of the curvature of the condensate) 



Rc= [Vo/n / d d rMr)[-Kr (r)]n(r) 



-1/2 



(4.6) 



d d r\^ \ 2 exp(ikr)-Mk) / d d r%{r)\Mr)\ 



and 



C(k) 

(4.7)' 

Note, that only requirement on the Hamiltonian of the 
system, that is needed to derive the inequality (|4.5|) is 
that it should not contain any explicit velocity depen- 
dence in the interaction and external potentials. 

Since in 2D scales as n, this case can be consid- 
ered as marginal and the condensate still can emerge 
even in an interacting system. This is because the usual 
log divergences inherent for 2D are cut off by a trap. 
The inequality (|4.5|> is a geometrical equivalent of the 
Bogoliubov-Hohenberg inequality, since it gives the lower 
bound for the ratio of the effective radius of the conden- 
sate to the de Broglie wavelength X^b ■ The second term 
of rhs (|4.5() can be used to obtain an upper limit on the 
possible condensate fraction as a function of temperature. 
Concrete examples of the application of (14.511 to_quasi-lD 
systems are given in the paper bv iFischerl ((2002) . 



One can also approach the problem of the condensate 
formation by directly analyzing the pha se fluctuations 
of the order parameter (for a review see iHellwee et all 
(1200 ll^ . Phase fluctuations are caused by the thermal 
excitations and are always present at finite temperatures. 
Note, that at very low temperatures density fluctuations 
in equilibrium are suppressed due to their energetic cost 
and can therefore be ignored. This assumption is not 
valid in the vicinity of a vortex core, but at very low 
temperatures, the vortex formation is negligible. 

As an aside, we mention that the concept of phase 
in quantum systems, introduced by Dirac as a canoni- 
cal conjugate observable to the number operator h, re- 
mains a controversial issue in certain circles. Formally 
it is known, that if h is an operator with a purely dis- 
crete spectrum (which is always true for the number 
operator), then there can exist no operator 9 such as 
the commutator [n, 9] = il holds. Different versions of 
phase-related operators have been constructed in order 
to overcome this difficulty (s ee for example the review by 
ICarruthers and Nietol (1968) and the t extbook on Quan- 
tum optics bv lMandeland" WolJ {19951) 1 . Alternatives to 
conventional symmetry breaking approache s have even 
been proposed (see the paper of IStenholml (|2002i1 and 
references therein ). An intriguing suggestion that the 
interference patterns of two atomic condensates can be 
explained with out ever evoking the notion of phase was 
put forward bv l Javanainen and Yol (Il99rl . 

In the present article we adopt the "conventional" and 
certainly more convenient approach, according to which 
the bosonic field operator takes on the form 



tp(r) = y/n (r) exp[i0(r)] 



(4.8) 



for the large number of particles. Here 9(r) is the opera- 
tor of the phase fluctuations and uq (r) is the condensate 
density at T = 0. 

To proceed with calculations it is convenient to ex- 
pand the phase operator in terms of the creation and 
annihilation opera tors for Bogoliubov quasiparticles (see 
iShevchenkol <ll992l) 1 



2Vn (r) 



k 



{(uk + v k )a k + (u k - v k )al) , 



where a k is the annihilation operator for the Bogoli- 
ubov excitation with energy e k , and u k , v k are excita- 
tion functions, determined by a bosonic equivalent of the 
Bogoliubov-de Gennes equations (for a general reference 
see the book of lde Gennea ((1966) ). Expression l|4.9fl can 
be obtained in the formalism of Bogoliubov transforma- 
tion generalized to an inhomogeneous case. 

Phase fluctuations in a quasi 2D system can be an- 
alyzed within the formal ism of the one-particle de nsity 
matrix (see the works bv lPetrov etafi ()2000Ll200lM 



(V> + (r)V(0)) = ^o(r)M0)ex P H(A#(r)) 2 )/2)] 



(4.10) 
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One should mention that the quasi two-dimensionality 
of the system implies that the scattering of particles ac- 
quires a 3D character, while the kinetic properties of the 
gas remain two-dimensional. 

It is clear from (I4.9|l that the estimation of the phase 
fluctuations (A6*(r) 2 ) requires a knowledge of the Bogoli- 
ubov quasiparticle s pectrum i n the i nhom ogeneous sys- 
tems (see papers of IStringaril lll996[) and lOnbe rg et all 
lll997l) and references in papers bv iPetrov et al\ (2000, 
2001)). This spectrum is discrete for T <C [i and for 
r>/j one can use the local density approximation. In 
the Thomas-Fermi regime for T \l one obtains the 
following approximation 



slow modes 



fast modes 



(A6(v) 2 ) ~ T\n(R/\ dB ). 



(4.11) 



Note, that (|4.11() does not depend on precise expression 
for the repulsive coupling constant. 

From (|4. lift one can estimate the characteristic ra- 
dius Rg of phase fluctuations (the characteristic length 
at which phase changes by 27r) to be Rg ~ XdBexp(Tg/T) 
with ksTg = N(haj±) 2 j \i. We thus arrive at the con- 
clusion that at low temperatures T <C Tg the character- 
istic radius of the phase fluctuations is larger than the 
size of the trap Rg 3> R±, so a true condensate exists. 
The emergence of a true condensate is attributed to the 
weakening of phase fluctuations induced by a trap, which 
introduces a low momenta cut off into the excitations in 
the system. At higher temperatures, T 3> Tg the system 
is characterized as a quasicondensate (Rg -C R±). 

The crucial effect of a trap f or 2D Bose gases was 
also emphasized bv iHo and Mai l)l999ft . They pointed 
out that long wavelength quantum fluctuations will be 
partially supp ressed due to th e gapped spectrum of col- 
lective modes l|Stringarill996j) and off-diagonal order will 
survive in 2D. 

Since there is an experimental evidence in support of 
BEC existence in 2D, the discussion is not yet closed. 
Quantum Monte Carlo simulations for bosons in a two- 
dimensional harmonic trap do indeed show that a signif- 
icant fraction of the p articles is still present in the lowest 
state at low energies (jHeinrichs and MullinLll99^ . 



C. Popov's approach 



In this section we consider the wav lPopovl l)l983j) gen- 
eralized the field-theoretical methods, developed by Beli- 
aev to finite temperatures. It is curious that the method, 
suggested by Popov in 1965, is conceptually very similar 
to renormalization-group approach, successfully applied 
in the 1970th to phenomena, un accessible to pe rturbative 
methods, such as Kondo effect teewsodll993|) . 

As usual one starts with the introduction of the tem- 
perature Green's function 



G(x,t;x,t) = — < ip(x, T)ip(x', t') > 

J e s ^p(x, r)tp(x', r')d^dip 

j e s dipdip 
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FIG. 5 Energy scales in Popov's approach: T/ln|eo/V| <C 



k'o/2m ~ T < (k' ) 2 /2m < e = 1/ma 2 



where S is the classical action of the Bose gas 

r 8 r _ r 13 

S= dr d 3 xip(x,r)d T ij(x,T) - / cLtH(t) (4.13) 
Jo J Jo 

and 



H(t)= J d 3 xiP(x,T)(-— -h)$(x,t)+ 

(4.14) 

Next step is the construction of the perturbation the- 
ory and corresponding diagrams arising from integrals 
of the type (|4.12l) . by performing the usual "trick" of 
separating out the condensate operators (|3.9|) . However, 
in case of Bose system the perturbation series converges 
very poorly for small momenta and frequencies. In other 
words, the infrared asymptote of the Green's function is 
singular. In order to avoid these difficulties, Popov sug- 
gested following modifications: the bosonic fields tp 



ip(x,r) 



t y exp[— i(kx — u!T)]a(k, ui) (4.15) 



(4.12) 



is divided into a short wavelength "fast" component ipi 
and a long wavelength "slow" component ^0 ('<P( x t t ) = 
ipo(x, r)+ipt(x, t)) (see Fig. 5). The momentum ko which 
separates the slow modes from the rapidly oscillating 
modes depends on the particular Bose system and only 
its order of magnitude can be estimated. Introduction of 
ko removes the divergences at small momenta, regulariz- 
ing the perturbation theory. 

A method of successive integration, first over the 
"rapid" and then over the "slow" fields, is then applied, 
using different schemes of perturbation th eory at differe nt 
stages of the integration (see chapter 4 in lPopovl |)1983|1 ). 
The fast modes "see" the slow modes as an effective con- 
densate ("bare" condensate according to Popov) with a 
superfluid density po = \ipo\ 2 - In Appendix II we give a 
succinct derivation of main Popov's results. 

This method of subsequent integration, developed by 
Popov, allows to estimate the low temperature asymp- 
totic behaviour of the one -particle Green's function, and 
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to derive a power-law decay of G(x,y) ~ \x — y\~ a for 
|x — y\ — ► oo (in ID and 2D) rather than the exponential 
decay that occurs at high temperatures. In 2D, as we 
have mentioned in HV.AI this signals the development of 
topological LRO at low temperatures. 

The analysis is based on the t-matrix description of 
the effective interactions, and the key property of the 2D 
t-matrix is that at low energies it vanishes, and at high 
energy cut-off the t-matrix diverges (see Appendix II). 
This results in an extremely small critical temperature 

4 In m(e //u) 

where e is a high energy cut-off and p is a chemi- 
cal potential. Bear in mind, that this is a mean-field 
derivation and the condition for the superfluid transi- 
ti on was assumed to be p = p n , because Popov (as well 
as iBerezinslril (|l970tll97i|) ') thought that at the critical 
temperature T Cl the superfluid density vanishes. 

The applicability of Popov's mean-field description is 
based on the assumption of a very small exponent a. For 
large a the probability of creation of quantum vortices 
becomes big and even this modified perturbation theory 
is invalid (see also discussion in llV.DI and the "corrected" 
many-body mean-field theory m llV.Efl . The applicability 
of the Popov's t-matrix description and the diluteness 
condition, derived by Fisher and Hohenberg is the main 
subject of Section llV.DI 



D. Diluteness condition and validity of t-matrix 
approximation 

We have already discussed that the perturbative treat- 
ment of the dilute weakly interacting Bose gas amounts 
to replacing the real potential by an effective two-particle 
t-matrix, obtained by summing up all ladder diagrams. 
From this point of view the diluteness condition deter- 
mines the range of validity of the t-matrix approxima- 
tion. 

An explicit form for the diluteness condition of 2D in- 
teracting Bose gas at finite temp e rature s was first intro- 
duced bv iFisher and Hohenberd l)l988lh They pointed 
out that singularities inherent to 2D systems (vanishing 
of scattering t matrix at zero temperature and classical 
divergence of phase fluctuations) might lead to drastic 
modifications of the usual dilute gas expansion. 

As we have discussed (see section lHI.B|) at zero temper- 
ature in 2D the diluteness condition na 2 <C 1 is replaced 
by 

nrf^ « L ( 4 - 17 ) 

ln(l/na ) 

Popov's theory can be used to demonstrate that 
at finite temperature, the above condition (|4.17|) 
is replaced by an even m ore stringent inequality 
ijFisher and Hohenberd. Il988ft 



Fisher and Hohenberg (FH) provided a heuristic 
derivation of this result, based on the Bogoliubov quasi- 
particle picture. Their analysis is based on the simple 
observation that the usual Landau quasiparticlc formula 
for the superfluid density 



p 



= l - 



d d k 



Kk 



pdj (2tt) 2 ' v (e«* - l) 2 



(4.19) 



where d is the dimension, does not have any singularities 
for d = 2, except in case when the chemical potential 
is small (p is introduced in (|4.19|l via the Bogoliubov 
quasi-particle spectrum £ 2 = n^l/ok 2 /m + A: 4 /4m 2 = 
pk 2 /m + k 4 /4m 2 . The validity of th is approximation 
for p is discussed in l|Beliaevl Il958albf )h By introducing 
the infrared cut off (ko ~ ^J~p) via the ansatz 



/'■ 1 



\ln(a 2 p) \ ' 



(4.20) 



the regularization of the integral of Ea. H4.19|) can be 
achieved, and one arrives at Popov's equations for the 
"superfluid" and "normal" densities 1|5.14[1 - 1|5.15[1 . 

The analysis of the temperature dependence of the su- 
perfluid density allows one to separate out three char- 
acteristic regimes: (i) the low temperature region, the 
physics of which is defined by phononic behaviour of the 
quasiparticles, leading to a superfluid density which de- 
pends on temperature as (1 — aT 3 ); (ii) a free parti- 
cle region, where p s behaves linearly with temperature; 
(iii) and a critical region, determined by the fluctuations 
around the critical temperature T c ( p s vanishes at T c ) 
(see Fig. 6.) 
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FIG. 6 Schematic phase diagram of a uniform dilute Bose 
gas. p a is the superfluid density, normalized by the mass 
density of the gas p = mn, Tkt is the critical temperature 
of the Kosterlitz-Thouless transition, Ta/f is the mean-field 
temperature, calculated perturbatively. Size of the transition 
critical region is defined by a parameter 1/ lnln(l/na 2 ), where 
a is characteristic s-wave scattering length. The region, dom- 
inated by phononic quasi-particle behaviour is of the width 
l/ln(l/na 2 ). 

The diluteness criterion is determined by the condi- 
tion that the critical region is small enough so that all 
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three regimes can be well separated. The width of the 
first regime is in fact given by Schick's small parame- 
ter l|4.17|l . while the size of the critical region is char- 
acterized by the double log l|4.18|l . The problem how- 
ever, is that for all practically relevant situations, even 
for very small na 2 , Fisher and Hohcnberg's small param- 
eter 1 / lnln(l / na 2 ) is still orders of magnitude greater 
than 1/lnil/na 2 ). This means that in practice the criti- 
cal region associated with Kosterlitz-Thouless transition 
is so large, that mean-field based approaches do not give 
any reliable results. Note, that the double log result was 
also reproduced by Fisher and Hohenberg in a more accu- 
rate way within a renormalization group treatment of the 
same problem. They have also estimated the superfluid 
transition temperature, which reads 



c mlnln[l/(na s ) 2 ] 

The results of FH work ijFisher and Hohenberst [19881 
have been confirmed in other approaches, see for exa m- 
ple virial expa nsion of a dilute Bose gas bvlRenl ll2004l) or 
RG analysis bylKolomeiskv and Stralev l lll992alrJ) and by 
ICrisan et all I^OOlj) . iPieri et al\ (|2001albT) demonstrated 
by analyzing the normal state by standard diagram tech- 
nique that the transition temperature 14.21|) appears as 
a lower bound for the validity of the i-matrix as a con- 
trolled approximation for the dilute Bose gas. 

The FH diluteness condition l|4.18(l is extremely strin- 
gent, and if straigh tforwardly applied to experimentally 
releva nt situations ijGorlitz et all l2001i iRvchtarik et all 
would mean that the systems observed to undergo 
a BEC phase transition in 2D are not actually dilute, 
and could nev er be so. This line o f reasoning motivated 
Liu and Wen llLiu and Wenl . l2002t) to come up with an 
exotic alternative scenario involving a two-dimensional 
strongly-correlated spin liquid. 

The extreme conclusions drawn from the FH diluteness 
criterion are nevertheless related to the general draw- 
backs of the Popov approximation. We will see in the 
next section, that in a more realistic model, which takes 
into account interactions in a self-consistent way, the di- 
luteness condition becomes much weaker. Moreover, in 
view of our previous discussion about the inapplicability 
of arguments based on homogeneous systems in the ther- 
modynamic limit to trapped gases, it would seem that 
the FH diluteness requirement is not really relevant for 
the experimental situation of the Bose gas in a magnetic 
trap. 

E. Other approaches: RPA, Many-body t-matrix, Monte 
Carlo 

In this section we review a range of diagrammatic ap- 
proaches that have built upon the early RPA and t- 
matrix approximation in order to improve the description 
of the superfluid of BKT transition and also the numer- 
ical methods, which allow to directly probe the critical 
region of the 2D transition. 



The first finite temperature generalization of Bo- 
goliubov rando m phase app r oxim ation i|lll.A[) was 
introduced b y jTserkovnikovl |)1964|) (english version 
iTserkovnikOTl l)l965j) ). His concern was to calculate 
the finite-temperature correction to the condensate den- 
sity in 3D dilute Bose gases with weak interactions. 
Tserkovnikov assumed that the average single particle 
kinetic energy is small, compared to the potential energy 
for all temperatures below T^ D . He also remarks that 
his approximation does not meet the Landau superfluid- 
ity criterion and that more precise equations should be 
sought in future work. 

The RPA method was furth e r dev eloped in the pa- 
pers of ISzepfalusv and Kondorl l)l974|) . whose main in- 
terest was the investigation of dynamics of the second 
order phase transition. Around the sam e time a lar ge- N 
approach was app lied to the Bose gas bv lAbd l)l974|) and 
lAbe" and Hikamil l)l974|) . who calculated the dynamical 
scaling for one-particle Green's function up to 0(1/N). 
Here the idea of the large- N approach is to expand the 
number of independent components of the Bose field from 
unity to N/2 using l/N as an expansion parameter. 
To produce a controlled large N limit, the interaction 
strength is scaled to be of order l/N. 

The RPA large-A^ method has been applied to 2D 
Bose gas bv lNoeueira and Kleinertl l)2006j) . The interac- 
tion in their approach is approximated by a 2D coupling 
constant, derived in the i-matrix approximation, consid- 
ered by Popov and Schick (see sections IIII.BI IIV.GI and 
lIV.Dfl . It is however known that in a large- N approach 
one can not simultaneously account for both particle- 
hole channel (RPA) and particle-particle channel in a 
well-controlled fashion. Nevertheless, the authors claim 
that the diluteness condition leads to the appearing of t- 
matrix diagrams first, while the n ext class of diagrams are 
those from particle-hole channel ijNoeueira and Kleinertl 
2006). This approximation results in the Bogoliubov 
quasiparticle dispersion containing a log correction due to 

low dimensionality & = ^ e 2 k + 2g 2 Dne k [l - ^ln(ka)], 

so that excitations in the system exhibit a roton-like min- 
imum. Note that the excitation spectrum is calculated 
under the assumption that one-particle Green's function 
and density correl ators share the same pol es (this prop- 
erty was derived bv lHohenberg and Martini ljl965() in case 
of 3D condensed Bose system). Would be interesting to 
check if these RPA results are confirmed in other ap- 
proaches. 

We now proceed to a discussion of the various gener- 
alizations of the two-body i-matrix approach. Though 
simple and elegant, the perturbative two-body £-matrix 
approach does have its drawbacks. The main problem 
is related to its inability to properly describe the criti- 
cal region in low dimensions. For example, the t-matrix 
method does not predict the Nelson-Kosterlitz universal 
jump in the superfluid density. In 3D the two-body t- 
matrix approach leads to a first order phase transition for 
the condensate density, which is the consequence of non- 
self-consistency of this first order perturbative approx- 
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imation (see also iGriffinl <ll988ft : E ee and Yand Jl958); 
iReatto and Stralevl l)l969|) ). 

Many of these problems can be solved if the many- 
body corrections, arising due to the surrounding gaseous 
medium are taken into account. This is the key idea 
in the "many-body t-mat rix approximat i on" (s ee the 
comprehensive review by IShi and Gr iffin (1998), low- 
dimensional systems within many-b ody t-matrix ap- 
proach are analyzed in the p a pers by | Al K hawai a et all 
l)2002j) : lAndersen et oZJ (|2002^ : IStoof and Biilsmal l)l993|) . 
for Hartree-Fock -Bogoliub ov study of a two-dim e nsiona l 
gas see recent works by iGies and Hutchinson! l)2004|k 
lOies et all 112005112004^ . 

Since the "many-body i-matrix" methods are exten- 
sively discussed in the literature, here we restrict ourself 
to a brief description of the method, providing all rel- 
evant references. The Bogoliub ov-Hartree-Fock (BHF) 
approximation (see lGriffinl l|l99(j) and the analysis of ex- 
citatio ns in a trapped 3D gas paper bv lHutchinson et all 
lll997M has a Heisenberg equation of motion for a Bose 
field operator of the kind l|3.24|l as a starting point 



, d , 
di ^ ' 



2m 



V +V ext (r) -// ^(r,i) + 



+ gip + (r,t)i(>(r,t)iP(r,t). 



(4.22) 



A short-range interaction is assumed among the atoms 
U(r — r') = gS(r — r'). Treating the interaction term in 
this equation in the self-consistent mean-field approxima- 
tion, one arrives at the equation 



2m 



V 2 +V e xt{r) -n 



+g{n (r) + 2n'(r))0(r) + gm' \r)<f>* (r) = (4,23) 

where no is the condensate density, n'(r) — (r)ip' (r)) 
and m'(r) — (?p' (r)ip' (r)} (anomalous average). In order 
to describe excitations in the system one should also write 
down the equation of motion for -0' 

+2gn(r)ip'(r,t) + gm(r)tj/ + (r,t) (4.24) 

(with n(r) = n + n' and m(r) = (f> 2 + ml). Equations 
(|4.23l) and (|4.24|) correspond to the Bogoliubov-Hartree- 
Fock approximation: a bosonic analog of the finite- 
temperature Bogoliubov-de Gennes equations. These 
equations can also be re-expressed in terms of a Green's 
function formalism. Note, that the appearance of anoma- 
lous averages in BHF formalism (which are not present 
in Popov approach), leads to a gap in the quasiparticle 
excitation spectrum. 

Many body effects are also effectively treated in 
a variational method , app lied to dilute Bose gases 
by iBiilsma and Stool l)l997[) . and many-body t- matrix 
methods whic h have been significantly improved in re- 
cent years (see lAl Khawaia et all l)2002|) ; [Andersen et al\ 



(2002)). A time depende nt BHF approximation has re- 
cently been developed bvlProukakis et all (|l998j) . Here, 
the authors l)Proukakis et adll998|) claim that the pseu- 
dopotential approximation U(r — r') = gS(r — r') should 
be imposed only after the effective interaction is ex- 
pressed in terms of many-body i-matrix. Both ap- 
proaches, the one used by Griffin et al. and another, 
developed by Stoof and collaborators, are qualitatively 
similar, in that they treat interactions in many-body t- 
matrix approach, but they differ in some details, for in- 
stance in selecting out the important diagrams. 

Let us now discuss some of the results of the many- 
body t-ma trix method for 2D syst ems. More than a 
decade ago IStoof and Biilsmal l)l993|) demonstrated that 
the infrared divergences, appearing in the two-body t- 
matrix treatment, can be elegantly eliminated when the 
effects of surrounding gas are taken into account. With 
this approach the universal jump in supcrfluid density, 
predicted by Nelson and Kosterlitz can be reproduced. A 
weaker diluteness condition, namely that of Schick l|4.17|l . 
defines the applicability of this many-body i-matrix ap- 
proximation, one that is satisfied experimentally in sys- 
tems, such as s pin polarized atomic hydr ogen, absorbed 
on 4 He surface l)Stoof and Biilsmal fl993) . 

The conclusions o f lPetrov et all (|2000ll200l have been 
confir med in recent investigations bv | Ande rsen et oT 



2002) and IGies and HutchinsorJ (I2004D: IGies et al. 



(2( 
(2( 



200a l2004j) . The theory of lAndersen et all l|2002l) . free 
of infrared divergences in all dimensions, allows for calcu- 
lation of the density profile of a (quasi)-condensate cloud 
of a gas for any aspect ratio of the trap (within local 
density and Thomas-Fermi approximation). At very low 
temperature, depending on the trapping geometry, the 
presence of a true condensate in t he equilibrium state is 
found. Hutchinson and cowo rkers ijGies and HutchinsorJ . 
12004 IGies ei a/1 12005L 12004^ also "see" within their HFB 
approach a macroscopic occupation of the ground state 
at low temperatures, implying the presence of a conden- 
sate state. 

To conclude, the presence of the trap appears to stabi- 
lize the condensate against the long-wavelength fluctua- 
tions and the BEC state can form at finite, though very 
low temperatures, when the discrete nature of the energy 
spectrum is taken into account. 

Most reliable description of the 2D Bos e gas to date 
is pr o vided by Monte Carlo simulations llKaean et al. , 
2000[ iProkof'ev etaH l200lt iProkof'ev and Svistunov . 
2002), because it allows to study the critical region of 



the BKT transition, which is effectively very large and 
therefore unaccessible to perturbative methods. The nu- 
merical analysis is simplified, by the fact that the critical 
properties of all XY models are the same (see HV.Bll . It 
suffices therefore to study the classical l^l 4 model on the 
lattice within a Monte Carlo algorithm. 

Consider, for instance, the temperature-dependence 
of the particle density in the critical region of the 
BKT transition, which follows from perturbative analysis 
l|Fisher and HohenbereLll988HKagan et q/.lll987llPopovt 
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Il983j) of weakly-interacting system 



where 



mT ln G 

2t: mU e f f ' 



(4.25) 



where U e ff is an effective interaction, proportional to 
/q D (|3.19|) . and C is the constant, which is not possible 
to eva luate within perturbat ive expansion in powers of 
U e ff l|Prokof'ev et all l200ll) . Monte Carlo estimation 
gives C — 380 ±3; this large value of C makes it virtually 
impossible to reach the limit of small U e ff for weakly- 
interacting system. 

At transition we obtain an an accurate microscopic 
expression for the critical temperature of BKT transitio n 
l)Prokof'ev et all l200ll IProkof'ev and Svistunovl |2002|) . 



T 



KT 



2tt 



n mhx(C/mU e ff) 



It is interesting to compare this density n to the quasi- 
condensate density n q and superfluid density n s in the 
critical region. It turns out, that n q /n is of order of 
unity, unless mU e ff is exponentially small, while the ra- 
tio n q /n s is of order of 2, which means that superfluid 
density is substantially smaller than quasi-condensate 
density at the transition. 

The temperature behaviour of various densities, ob- 
tained in Monte Carlo procedure, can be used for check- 
ing whether RG and perturbative approaches essentially 
overlap. Indeed, Monte Carlo simulations have been 
able to capture the crossover between the mean-field 
behaviour and the crit ical fluctuation region describe d 
by the KT transition Jprokof 'ev and Svistunovl |2002|) . 
IProkof'ev and Svistunovl l)2002|) show that this crossover 
is characterized by a universal ratio of the superfluid and 
quasi-condensate density. One can also see that the con- 
ventional mean-filed result n s /n ~ 1 — T/Tkt is not 
valid anywh ere, while the modified mean-fie ld theory in- 
troduced bv IProkof'ev and Svistunovl l)2002|) can predict 
accurately the behaviour of the quasi-condensate density 
up to Tkt- 



F. Breathing modes of 2D systems 

At the end of this chapter we consider a universal prop- 
erty of a two-dimensional gas with a contact interaction, 
confin ed in a harmonic potential. iPitaevskii and Roscbl 
i|l997|) predict that such a system develops oscillations or 
breathing modes, which can be probed experimentally or 
in simulations and thus can serve as a practical criterion 
of the two-dimensional nature of a system. 

The appearance of breathing modes is related to 
a hidden "Lorentz" symmetry inherent to any two- 
dimensional Hamiltonian of the following general form 



i v 7 i<j 



(4.28) 



and H ex t = ^muu 2 r 2 is a harmonic potential. 

It is readily seen that Hq is scale invariant in case of a 
local 2D interaction 



U(Ti 



9 



(4.29) 



H = H + H e 



(in fact it is scale invariant for any potential with a prop- 
erty U(lr) — U(r)/l 2 ). The presence of a trap breaks the 
scale invariance of Hq. Note, that in principle the scal- 
ing invariance of the Hamiltonian Hq is broken in 2D, 
because then scattering phase shift is energy dependent 
due to the logarithmic dependences characteristic of two 
dimensions (phase shift is proportional to the coupling 
(4.26) constant g 2D or to \/ln(ka)). The energy-dependent 
phase shift sign als the breaking of scale invariance at the 
quantum level ijCabo et all . Il998|) . But this symmetry 
breaking is explicit and is not attributed to any phase 
transition physics. 

In spite of the breaking of scale invariance, because of 
a special property of the harmonic oscillator, a power- 
ful spectrum generating symmetry still exists. That can 
be seen from the commutator [H ext ,H] — ioJ^Q, where 
Q = l/2X)i(Pi r i + r iPi) is the generator of scale trans- 
formations. One can check that [Q, Hq] = 2iHq and 
[Q, H ext ] — —2iH ext . These results can be formulated 
within the well-known algebra of SU(1,1) or SO(2,l) sym- 
metry group, i.e. the 2D Lorentz group. 

Starting from the lowest energy state eo one can pro- 
duce higher order states with energies eo + 2nu>o (n = 
1, 2, ...) by applying one of the SO(2,l) group generators 
L+ = (Li + iL 2 )/V2 where L\ = (Hq - H pot )/2uj Q , L 2 = 
Q/2 (the correspondent annihilation operator is L~ = 
[L\ — iL-i)j\pi). The excitations with the energies 2ruvo 
are associated with the breathing, or pulsating modes of 
the system. 

As an example, the authors ijPitaevskii and Roscbl 
considered the classical Gross-Pitaevskii equa- 
tion and predicted the existence of undamped breathing 
modes in the condensate. The appearance of transverse 
breathing modes with a frequency equal to an integer 
multiple of the trap oscillation frequency was observed 
experimentally in an elongated condensate of 87 i?6 atoms 
l|Chevv et all 12002]) . Numerical simulations (exact diag- 
onalization) seem to indicate the existence of dipole or 
breathing modes i n a 2D system even for relative ly small 
number of atoms l|Haugset and Haueerudl il998|) . 

BHF study of a 2D Bose gas by G ies et al. 

(|Gies and Hutchinsonl I200I iGies et a71 120051 1200^1 has 
also showed that at low temperatures the frequency of 
the lowest lying excitation (n = mode) is precisely 2wo, 
independent of the interaction strength. At high tem- 
peratures the frequency of this mode shifts to a lower 
frequency region, being modified by the addition of the 
potential from the static thermal cloud. 



(4.27) 
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Nevertheless, it is important to note that a (5-function 
is not well-define d in two dimensions due to l ogarithmic 
UV divergences l|Pitaevskii and R.oschl Ii997), that are 
cut off by the finite range of interaction, whether it is a 
small or large effect, should be investigated. 



V. CONCLUSIONS AND OPEN QUESTIONS 

We have surveyed a number of theoretical issues arising 
in the field of a weakly-interacting uniform or confined 
in a trap dilute Bose system at low temperatures in 2D. 
The underlying physics of such a system depends on the 
size of the system, degree of its inhomogeneity, and tem- 
perature. 

If the system is very large and uniform, one might ex- 
pect realization of the BKT transition, characterized by 
the presence of a topological order below the critical tem- 
perature Tkt down to zero temperature when the true 
long-range order (BEC) forms. Perturbative approaches, 
based on a low density approximation and point-like or 
short-range interactions, surveyed in this Colloquium, 
are not really suited to a description of the vortex exci- 
tations of the ordered phase of BKT transition, however, 
these methods provide a good description of many phys- 
ical properties. For example, a modified version of the 
mean-field theory, the "many-body i-matrix approach" is 
able to capture the second order nature of the phase tran- 
sition in 3D and the Nelson-Kosterlitz universal jump of 
superfluidity in 2D. 

In experiment the low temperature regime of a quan- 
tum gas is achieved by confining the atomic system in 
an external potential. The proper description of a prac- 
tically realized system requires therefore the inclusion of 
the inhomogeneities, introduced by a trap. In 3D the 
effect of the trap is not very pronounced, and it is possi- 
ble to calculate the correction to t he critical temperature 
due t o interactions perturbatively ijArnold and Tomasikl 
l200l|) . while in the un i form case the RG study is re- 
qu ired llLedowski et all \2004i ffo r details see the book 
bv lPitaevskii and Stringaril l)2003j) l. 

In 2D, as we have seen in the case of a gas without 
interactions, the presence of the trap dramatically mod- 
ifies properties of the system (density of states), so that 
BEC becomes possible at finite temperature. The in- 
clusion of interactions into the picture is a complicated 
task. First of all, the system is inhomogeneous and all 
previously developed perturbative methods, such as the 
t-matrix approach are strictly speaking not suitable for 
its analysis. In principle, one should solve the many-body 
scattering problem in a trap, taking into account the dis- 
crete spectrum and the finite range of a potential, which 
is extremely difficult. One can of course consider a sim- 
plified problem of a quasi-homogeneous trap and adapt 
well-studied techniques for that case. As we have seen, 
the main effects of the trap are indeed captured at least 
qualitatively within such a scheme. 

Intuitively, it is clear, that inhomogeneities would tend 



to suppress the universal jump of superfluidity, and would 
rather favor the "true" BEC state at low temperatures 
in a system with finite number of particles (because the 
long-range thermal fluctuations will be quenched by a 
trap). 

Very recently IHolzmann et all l|2005|) have analyzed 
the behaviour of the weakly-interacting trapped system 
in the thermodynamic limit within local density approxi- 
mation (LDA). They have shown, that although the uni- 
versal Nelson-Kosterlitz jump is indeed not present, the 
system does undergo BKT transition at the temperature, 
somewhat lower t han T„ 2D 112.1211. 

In a paper by ISimula et a l. (2005) the authors have 
predicted that both a BE condensed phase and a KT su- 
perfluid phase, separated by a first order transition, will 
be present in a 2D trapped gas. The authors arrive at 
their conclusions by comparing the Helmholtz free en- 
ergy of the ground state, characterized by a condensate 
wave-function and an excited state, containing a vortex 
pair in the ordering field. Based on entropy consider- 
ations, they find a critical temperature T c above which 
a thermally activated transition to the state containing 
vortex-pair excitations becomes favourable. These con- 
clusions require, however, solid confirmations from both 
theory and experiment. 

A challenging problem is however the actual estima- 
tion of the BEC transition temperature in a 2D interact- 
ing system and the relationship of the quasi-condensate 
state to the KT vortex-pair plasma state, depending on 
the geometry of the system and number of particles. It is 
also important to develop numerical methods adapted for 
experiment. One of the pr omising deve l opments in thi s 
direction has been made bv lDavis et aO ll2001albL 12002). 
who have proposed a projected Gross-Pitaevskii equa- 
tion formalism, which allows the efficient investigation 
of finite temperature properties of the equilibrium con- 
densate state - even in the region where the Bogoliubov 
theory fails. An extension of this method to harmoni- 
cally trapped c ondensates has been con sidered in the re- 
cent paper bv jB lakie an d David l)2005|) . Just this year 
ISimula and Blakia l)2006j) have analyzed the 2D Bose 
gas by classical field methods, adopted from quantum 
optics. They have demonstrated that two distinct su- 
perfluid phases, separated by thermal vortex-antivortex 
pair creation , exist in experimentally producible quasi- 
2D Bose gas. ISimula and Blakiel l|2006j) provided a strong 
evidence that a strange ("zipper") i nterference pattern 
observed in a recent experiment by IStock et al\ II2005D 
can be explained by a presence of a vortex excitations in 
an experimental system. 

One can thus tentatively identify the following phases 
in a trapped 2D Bose gas (Fig. 7): 

• phase I: low -temperature true BEC phase; 

• phase II: KT vortex-antivortex pair superfluid, or 
quasi-condensate, or condensate with a fluctuating 
phase, at the transition Tkt superfluid universal 
jump is almost suppressed; 
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FIG. 7 Schematic phase diagram of a 2D trapped weakly- 
interacting dilute Bose gas. po stands for the density of the 
condensate, p a is the superfluid density, pmf is the mean-field 
density, which can be estimated perturbatively. Tbec is the 
crossover temperature from the superfluid regime to the true 
Bose-Einstein condensation phase. Tbkt is the critical tem- 
perature of the Kosterlitz-Thouless transition. Tmf marks 
the critical region of BKT transition. 

• phase III: critical region of the BKT transitiona; 
vortex pair dissolve, above Tkt vortices are un- 
bound and free; 

• phase IV: above the mean-field temperature it is a 
normal fluid, no local order parameter or vortices 
exist. 

The phase diagram of a real system will depend on 
many factors, as we stressed at the beginning of this 
Chapter. Many aspects of the diagram, depicted on Fig. 
7 require careful investigation, and reliable confirmation 
from both theory and experiment. 

In ending we attempt to summarize some of the open 
questions: 

• Nature of superfluid phases of a 2D weakly- 
interacting Bose gas: what is the nature of the 
crossover to a superfluid phase? What is the ex- 
plicit relation between superfluidity and the "quasi- 
condensate" state? Is there a crossover to BEC 
state at low temperatures? and if yes, under which 
conditions? 

• Can we help experimentalists to "see" the vortex 
excitations in the superfluid state, to really identify 
the KT state? There is a clear need in good vortex 
detection methods. Can, for example, disorder help 
us to pin the vortices? 

• What are in general measurable physical prop- 
erties, which delineate between the coherent 
condensed state and superfluidity? Initial 
progress in this d ir ection has been already made, 
iPolkovnikov et all (12005^ have suggested how to 
identify the KT transition from experimentally 
measured interference pattern. 

• Can we solve the many-body scattering problem in 
the trap? Does diluteness of the gas simplify this 
problem? 



• Is local density a good approximation for descrip- 
tion of the experimental systems, and if yes, under 
which conditions? 

• Could one justify a large N approach which im- 
proves on the existing method by incorporating the 
i-matrix approximation? 

We h ave also discussed the dilute ness condition, de- 
rived bv lFisher and Hohenber d |l988) under certain con- 
ditions of the transition to superfluid state. In a finite- 
sized experimental system this condition is not really ap- 
plicable, and one should use the criterion of Schick l|4.17|l . 
which can be seen from the analysis of qua ntum fluctua- 
tions of the 2D BEC at zero temperature ijPetrov et all 
I2004D . 

We have considered only low density approximations. 
When the gas is dense, approaches such as that of Gross 
and Pitaevskii are not applicable. In such high density 
regimes, new "strong coupling" approaches are required. 
One of the possible solutions may be the slave-boson ap- 
proxima tion, which is valid for ha rd-core bosons at any 
density l)Zierier and Shuklal fl997^l . 

We have not discussed in our Colloquium the role 
of disorder in the continuum Bose system, though it 
could be a subject of a separate review and opens up 
a lot of interesting perspectives. Recent Monte Carlo 
studies predict, for example, that for strong disorder 
the system enters an unusual regime, where the su- 
perfl uid fraction is smaller tha n the condensate frac- 
tion ijAstrakharchik et all 12002). W eak disorder c an be 
treated within Bogoliub ov theory ijGiorgini et all 1 19941 
iHuang and MeneL Il992|) and the striking result of this 
is that disorder is more active in reducing superfluidity 
than in depleting the condensate. These results suggest 
that the relation of superfluidity and Bose-Einstein con- 
densation require further theoretical and experimental 
investigation. 

One can not but mention the rapidly increasing in- 
terest to the cold gases with dipole-dipole interactions, 
which are responsible for a variety of novel phenom- 
ena in ultracol d dipolar system s, see, for instance, 
iFischerl (120061): iPedri and Santosl (|2005(k ISantos et all 
l)2003tl2000|) : IStuhler et al\ l)2005l) and references therein. 

Finally, the problem of measurable quantities is in fact 
one of the most important in the context of trapped 
Bose gases. Unlike an electron system, one cannot attach 
"leads" to the system and measure transport properties 
of a condensate cloud. One of the most pressing practical 
needs for theorists and experimentalists is therefore the 
development of controllable new methods to probe the 
trapped condensate. 



Appendix I (for section HV.BI 

Sketch of the derivation of Mermin-Wagner- 
Hohenberg theorem 
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One should use the Bogoliubov inequality 



the commutation relation 



1 



({A,A + })([[C,H],C + }) >T\([C,A})\' 



(5.1) 



[<ty(r), <ty+(r')] = 6(r - r') - Mr)VoV), ( 5 - 5 ) 



where the average (X) = Tr (Xe"^ ff ) /Tre~@ H , and op- 
erators ^4 and C are such, that the ensemble averages on 
the lhs of (|5.1(l exist. The inequality (|5.1(l follows quite 
straightforwardly from the Schwartz inequality 



(A,A)(B,B)>\(A,B) 



(5.2) 



where a scalar 



product is defined by (A, B) = 
where xab is the Fourier transform of 



the response function XAB(rt,r't' ) = (■^■\A(rt) 1 B(r't')]) 
(see for example the textbook bv iForsterf <f 



1199(1 ). 

In the c ase o f a Bose system the derivation of 
iHohenberel l)l967|) is based on Bogoliubov and Schwartz 
inequalities and the / sum rule 



duil k 2 n 
T I XAA+(~kuj) = . 

IT U! m 



(5.3) 



Derivation by Fischer of geometrical analog of 
the Hohenberg inequality 

The bosonic field operator is as usual decomposed into 
condensate and noncondensate parts 



■0(r) = ip (r)a + 8ip(r). 



(5.4) 



The key observation of lFische'rl J2002) is that the Bogoli- 
ubov prescription should be applied after implementing 



for otherwise, the second term on the rhs of (|5.5|) . which 
turns out to be crucial for calculating the condensate 
fraction correctly, would vanish. 

The operators A and C in the Bogoliubov inequality 
(|5.1() are chosen to be smeared excitation and total den- 
sity operators 



d d rf A {v)5^{v)- C = d d rf c (r)Sp(r), (5.6) 



.4 



where f A and fc are carefully chosen "smearing func- 
tions" (/ c (r) ~ VoW, f A (r) ~ e* r ). Next, the /- 
sum rule analogous to l|5.3|l , can be derived in coordinate 
space. 



Appendix II (for the section HV. CI 

Popov's approach 

To derive the phase transition curve for a two- 
dimensional interacting Bose gas, one needs to explore 
the finite-temperature behaviour of the t- matrix l|3.14f) . 
The Bethe-Salpeter equation (Fig. 3) in the Matsubara 
representation reads 



r(pi,P25P3,P4) = C/ki-ka ~-gy^2 ^qG°(ki - q, iu)\ - lull) 

xG°(k 2 + q,iw 2 + iuJi)T(p 1 -q,p 2 + q;Pz,Pi), 



(5.7) 



where u)j = 2irjT is an even Matsubara frequency, and 
the four-dimensional vector pj = (kj , uij ) represents the 
momentum kj and frequency u>j of the particle before 
scattering (j = 1,2) and after (j = 3,4). Energy and 
momentum conservation requires pi + P2 = P3 + Pi- 

The main contribution to the sum over internal mo- 
menta in l|5.7f) comes from k ~ a~ 1 which is due to the 
potential, discussed above. Since a -1 ^> VT ^> .//I, 
the fi dep endence in the Green's function can be safely 
neglected (|Popovi Il983l) . Consequently, after integrat- 
ing over frequencies Eq. 1)5. 7JI is reduced to a £-matrix 
equation 

/dp* 1 
— — 2?7(pi -p')-pr, t(p',p 3 ,z) = 
(2tt) z p u /m — z 

= U( Pl - P3 ) (5.8) 
Schematically, the i-matrix equation can be expressed 



as t z + UR z t z = U, with R z = l/(p 2 /m - z). The 
operator (1 + UR Z ) can be inverted and we get U — 
t z {l — Rztz)^ 1 . In this fashion the interaction is elim- 
inated from the i-matrix equation and we arrive at the 
Hilbert identity t z —t Za = t z (R Z Q — R z )t ZQ . This last equa- 
tion is readily integrable, since at low energies (pi,P2 <C 
1/a, \z\ <C l/(ma 2 )) we can neglect the momentum de- 
pendence of the t matrix. The energy zq defines an ar- 
bitrary high-energy cut off of the order l/(ma 2 ), so that 
t(zo) 3> t(z), and we obtain the long- wavelength asymp- 
totic of the i-matrix in 2D 

47T 

mm[eo/i-z)] 

We see that in 2D, the i-matrix vanishes in the limit 
Pi j P2 j z —* and in fact diverges at the high-energy cut 
off e = \z \. 
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Next we need to integrate out the high-energy modes 
with momenta k > k' in our functional cxp S <|4.13[) . The 



cut-off k' is defined as 



max(\n\,T) « « e 

2 m 



(5.10) 



The result of this integration is the reduced action 



s , = 
1 



k 2 

. iu)-- \-fi ) a + (p)a(p) - 

£ — ' \ 2m 



2/3P 

Pl+P2=P3+P4 



t'a + ( Pl )a+(p 2 )a(p 3 ) a (p4) (5.11) 



where all the summations are cut-off at k — k' and the 
potential is replaced by a t-matrix with 

47T 



t' = t'{uj) 



TOln(e /[(fco 2 /"i) - iuj}) 



(5.12) 



Now the functional expS" is to be integrated over the 
variables a + {p),a{p) within the momentum shell fco < 
k < fcg, where kg is defined from the inequality fcg/m >> 
T/ln|eo/V| an d serves to distinguish between slow and 
rapid particles. Variables with momenta smaller than ko 
are taken into account in the action by the transformation 



a+(p),a(p) - (po(k Q )f3V) 1/2 6 pQ (5.13) 

here po{ko) is the density of slow particles and one can 
introduce the density of fast particles pi(k ). 

After the transformation (|5.13(l one can make use of 
standard perturbation theory formalism and derive ex- 
pressions for the densities po and p\. Their sum gives 
the total density p = po + p\ = p n + p s , which is inde- 
pendent of auxiliary momenta kg and k' Q 



mp n , - 1 



d 2 k 



k 2 



2me(k) 



n B {k) (5.14) 



where n B (k) — (e l3ek — 1) 1 , and the formulae for the 
normal and the superfluid component densities read 



Pn = 



J 3 

(2 
mp 

Ps = — (lne /M 

J2 



J d 2 k^-n B (k)(l + n B (k)) 

^(\ne /p - 1) - 
1 fd 2 k ]2 



--(3{n B {k) + l) 



. (5.15) 



Eqs H5.14|l and l|5.15|l define the thermodynamics of the 
system below the phase transition. Above the phase tran- 
sition one can use the standard perturbation theory for 
calculation of the Green's functions and the thermody- 
namical quantities. The critical temperature is derived 
now from the condition p = p n at the transition. 
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TABLE I Properties of the ideal gas. "DOS" stands for 
density of states p, "TDL" stands for thermodynamic limit. 
Tb ec is the critical temperature of Bose- Einstein condensa- 
tion depending on the dimension. JV is the number of par- 
ticles, V is the volume of the system, loho is the geometric 
average of the oscillator frequencies: in 2D cu ho = (uJxUjy) 1 ^ 2 
and in 3D ui ho = (tUxUjyLJz) 1 / 3 . 



ideal gas property 
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DOS 
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P2D ~ const 


P2D ~ e 
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TV — > oo, V — > oo 


3D: limw W) iV 1/3 = consf 
2D: limcj M 7V 1/2 = const 
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